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ABSTRACT 
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A missile system In which L laxinchers were available in 
a row and S spare missiles were stored in a continuation of 
this row was considered. Additionally, an assortment of 
warheads was available, but only one could be mounted in 
each missile. The ability of this system to meet a demand 
for a particular warhead was dependent on the mix of war- 
heads in the L+S missiles and the order in which the 
missiles with their varying warheads were mounted on the 
launchers and in storage. Since the missiles were in a 
continuous row, some missiles would not be available for 
firing before the missiles in front of them had been fired. 

Because of the myriad possible permutations of the L+S 
warheads, determination of optimal ordering of the warheads 
by simulation was found to be too time-consuming. Hence, 
analytic expressions were developed after an appropriate 
measure of effectiveness was devised. A computer program 
was developed to perform the recursive calculations 
necessary . 
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SYMBOLS 





number of targets missed due to 
availability between the ith and (i-l)st 
reload, in the section, in the k^h 

battery 


f f f 


number of targets missed after the storage 
barn was emptied, by the jth section, in 
the k^h battery 


AgjAe 


number of targets missed due to availability 
in a section and a battery respectively. 


B 


number of batteries in a defense 


®S»®B 


number of targets missed due to reloading 
in a section and a battery respectively 


PqCx) 


completely specified cumulative distribution 
function 


L 


number of launchers in a section 




number of targets missed during the i 
reload, in the jth section, in the k^h 
battery 


Ms,Mb,Mjj 


-total number of targets missed by a section, 
a battery, and a defense respectively 


^max 


maximum number of targets missed 


\in 


minimum number of targets missed 




probability of a request for a type j 
warhead 


^k’^kA 


weighting factor for a target missed due to 
availability by battery k, when it is of 
type A 


S 


number of missiles in storage in a section 


^k’^kA 


weighting factor for a target missed due 
to reloading by battery k, when it is of 
type A 


^15,000^^^ 


observed cumulative distribution function of 
a random sample of 15,000 observations 
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T 



number of reloads required in a section 



X., type of warhead stored on launcher i, 

^ i = 1,2,.. .,9. e {1,2, 3 , 4 , 5 } 

yi>y 2 >y -3 probability that a request is for the 

^ warhead on launcher #1, #2, or #3 

respectively. yi “ Pj 1^* = j , i = 1,2,3 
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I . INTRODUCTION 



Current articles on missile allocation problems make 
certain tacit assumptions. The primary assumption is that 
all missiles are mounted on their ovm launchers. Thus, the 
problem of determining an optimal mix of missiles becomes 
one of determining the optimal mix of all missiles emplaced. 
However, a question comes to mind: What if emplacing each 

missile on its own launcher is infeasible for reasons of 
cost, mobility, maintenance, etc.? Then the missile complex 
will consist of a certain number of missiles mounted on 
launchers with the remainder in storage . The problem then 
becomes two problems: Should the mix of missiles on 

launchers be the same as the mix of missiles in storage? 

If not, what should the respective mixes be? The two-part 
problem could be simplified back to the one-part problem by 
assuming that the time necessary to move the missiles in 
storage to launching positions will be so long as to render 
the stored missiles unavailable to any short run situation, 
that is, the storage is a reserve for future battles. 
However, if the system is set up so that the stored missiles 
can indeed be made available relatively quickly, the 
problem of respective mixes becomes crucial. 

A further complication also arises when the stored mis- 
siles' are not immediately available for reloading, but are 
placed in a sequential storage system. An analogy would be 
a magazine-fed rifle with a variety of types of rounds 
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available. The problem would consist not only of the mix 
of rounds In the magazine, but also the order in which the 
rounds are stored in the magazine. The order is critical 
because the rounds are fired sequentially. The situation 
can be understood by the following example: Suppose the 

magazine holds ten rounds and the only round of a certain 
type is placed tenth in the magazine. When the firer wants 
to fire that round, he must either fire or eject all of the 
intervening rounds. In either case, poor utilization of 
resources was caused by the ordering of rounds in the 
magazine . 

Returning to the missile problem, it is not likely that 
a variety of types of missiles would be stored or emplaced 
together. However, the same situation would arise if one 
type of missile could be fitted with one of a variety of 
warheads. Even though the missiles would have identical 
external configurations, they could be considered different 
missiles because their use would vary with the warhead 
carried. 

In summary, a missile site in which the missiles could 
be equipped with a variety of warheads and in which the 
missiles are placed in a ready storage area with a sequential 
reloading mechanism appears to present difficulties in 
allocation that are not present in systems commonly studied 
in the current literature. 

Consequently, this thesis will attempt to establish a 
methodology for evaluation of missile sites designed as 
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described above. The first step will be to establish a 
measure of effectiveness suitable for such a system. When 
this has been established, further questions can be 
answered. The first question- considered will be: Does 

ordering of warheads have any significant effect? If not, 
the question then becomes an academic exercise. If ordering 
is indeed significant (as it intuitively appears to be), 
further characteristics such as optimal ordering, optimal 
mix, or suboptimization, can be explored. 

Since a general system will be considered, no definitive, 
real-world answers will be obtained. However, by exploring 
different possibilities, directions for further research 
will be indicated. In short, this thesis will describe the 
necessity of a separate methodology for evaluation of 
sequential ordered systems. 
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II. THE PHYSICAL SYSTEM 



For this study, the launching area consists of three 
separate sections of three launchers and a storage barn 
holding six missiles. Thus the entire launching area (one 
battery) contains twenty-seven missiles when fully loaded. 
Further, it is assumed that five separate types of warheads 
can be installed in each missile [Headquarters, Department 
of the Army 1968], Additionally, the mix of warheads in 
this battery is considered to be determined by some other 
agency. Any warhead can be mounted on any missile with all 
missiles completely identical in external configuration and 
hence, no specific considerations need be given about placing 
a certain type of warhead in any particular one of the twenty- 
seven available missiles. However, once emplaced, each mis- 
sile has utility dependent upon which type of warhead is 
installed and where the missile is located in its section. 

An additional restriction upon availability of missiles 
is that a missile in a section can be moved to another sec- 
tion only at a time expense of several hours since the mis- 
sile has to be manually disassembled, moved to its new sec- 
tion, and reassembled. For safety reasons then, a section 
can not fire while a missile is being stocked. Another 
feature is that men must manually move the missiles from the 
storage locations onto the launchers during the reloading 
process. Hence, a section is out of action while either 
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reloading or restocking is accomplished. The launching area 
has enough personnel to conduct reloading operations in all 
three sections simultaneously. However, only one missile 
can be moved at a time in each section. The three sections 
are separated far enough that two sections can be reloading 
when the third section is firing. Since the restocking 
operation is so time consuming, its ramifications are con- 
sidered to be strategic, rather than tactical, and are not 
considered in this analysis. 

Because of the sequential organization, each missile 
has to be moved across all positions between their storage 
position and the destined launcher. This is best illus- 
trated by an example: 

IIHDCKK><><><>0 

123123^56 

I \ » — < 

Launchers Storage Barn 

Suppose the missile in storage #5 has to be moved to 
launcher H2. The missile has to move across storage posi- 
tions #4, #3, #2 and #1 as well as launcher #3- For this to 
be feasible, the five Intervening positions have to be empty. 
The only way this can occur is if the missiles had been 
fired, l.e. no means exists of disposing of a missile 
preventing a desired move. These missiles could have been 
fired only from the launchers. Hence, the four empty 
positions in the storage area indicated that previous to this 
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move, the first four missiles had been moved to the 
launchers and fired. 

Consequently, if a desired warhead is in storage #4 and 
storages and #1 are full, the desired warhead would 

be inaccessible and, in effect, nonexistent for the engage- 
ment. This inaccessibility is obvious when it is considered 
that at best (for this situation) all three launchers could 
be empty. If this was the case, the missiles in storages 
#3> #2, and #1 could be moved onto the launchers. But then 
the desired missile in storage still could not be fired 
because all launchers would be full. If one of these 
missiles was fired, the desired missile could then be moved 
onto a launcher. In summary, one missile has to be fired 
before the desired missile could possibly be fired. 

The decision as to which of the five types of warheads 
available is desired for each firing is made external to 
the launching area. Due to this, the decision in the 
launching area consists of either being able to satisfy the 
request for a warhead or not. This capability is dependent 
upon the total mix of warheads in the launching area, and 
the order in which the warheads (in their missiles) are 
placed in the three launching sections. It is also assumed 
that each request for a warhead is Independent of all others 
and that the Pj are available from the intelligence 
community . 

The over-all firing pattern for this system is a shoot- 
look-shoot technique. That is, after a missile is fired. 
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the results of Its detonation are assessed before the next 
round is fired. Thus the battery would have only one missile 
in the air at a given time. This has the effect that the 
battery engages only one target at a time no matter how many 
are within range. Thus, if many mobile targets are available, 
the battery would engage them serially until: (1) all targets 

are destroyed, (2) no targets are left within range, (3) no 
missiles are available on the launchers (reloading would 
enable reengagement), or (4) all missiles have been fired. 
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III. MEASURE OF EFFECTIVENESS 



A measure of effectiveness must concern the ability 
of a battery to perform its mission i.e. engage targets. 

This can possibly be measured by the number of targets not 
engaged out of the total targets available. This possibility 
has certain drawbacks. One is that the number of targets 
missed would be influenced as much by number of targets 
available as by the mix and ordering of warheads e.g., if 
fifty targets- were presented, at least twenty-three would be 
missed (only twenty-seven missiles available); if 100 tar- 
gets were presented, at least 73 targets or 73 ^ would be 
missed, etc. This difficulty is avoided by considering the 
number of targets not engaged prior to the last missile 
being fired. This end point is used because restocking is 
not considered possible in a tactical time frame. 

As was mentioned in the description of the system, only 
one target may be engaged at a time, regardless of the. number 
of targets available. Hence, if a sufficient number of 
targets are available, ordering of warheads and reloading 
time are the governing factors of number of targets not 
engaged because no time would be available for reloading 
due to lack of targets. Actually, this condition can exist 
for analysis by assuming that the opponent is using saturation 
tactics. This is not a really gross assumption; it is 
merely assuming that the opponent is using optimal tactics 
[Waddell I963]. This optimality is easily understood when 



one considers the situation that would exist if saturation 
tactics were not used. If the targets were spread out in 
time and/or space, the battery could possibly be able to 
reload without paying a penalty in terms of inability to 
engage targets due to reloading. 

A discussion of the role of time in this analysis is 
needed. When a section goes out of action in order to re- 
load, the length of time that elapses until the section is 
again ready to fire has a direct effect upon the engagement 
capability of the battery. In order to make the transfor- 
mation from distance each missile is moved during reload to 
time elapsed, it is assumed that the time required to move 
any missile to its adjacent position is constant; this 
constant is called one time unit . Hence to move a missile 
from storage #4 to launcher #2, five moves are required and 
consequently, five time units elapse. If, in addition, a 
move is made from storage #5 to launcher #3j five more time 
units would elapse or a total of ten time units would be 
required for the reload process. 

The next problem is to relate the reload time to tar- 
gets not engaged. As has already been explained, sufficient 
targets are available so that the battery can engage the 
next target immediately upon the detonation of the previous 
missile. Hence, the time between firings will be the fly- 
out time plus the reaction time, which is assessment time 
(was the target destroyed?) plus decision time (reengage 
the same target? what warhead?). A common assumption in air 
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defense analysis is that the reaction time is constant 
[United States Army Air Defense School 1971] and was assumed 
to be so in this paper. Consideration of the fly-out time 
immediately reveals that this Is a variable depending on 
speed of the missile, speed of the target, range of the tar- 
get, path of the target, etc. If this fly-out time plus 
reaction time for some firing Is longer than one time unit 
(as defined above), the battery could move one missile one 
place in each of the three sections without any decrease in 
engagement ability for it could not fire while the missile 
was in flight anyway. This would be an obvious benefit 
for the battery. However, virtually the only way to deter- 
mine various fly-out times during a battle would be simula- 
tion due to the myriad parameters mentioned above. In this 
paper, it is assumed that the fly-out time is constant in 
order to enhance analysis. 

Consideration of the magnitude of the constant fly-out 
time is necessary. Since both fly-out time and response time 
are constant, time between firings is now constant. If the 
time unit for reloading as defined above is less than this 
constant time between firing, the battery would be able to 
take advantage of the situation mentioned in the previous 
paragraph, i.e. movement of missiles with no penalty. 

In the situation in which the reload time unit is greater 
than or equal to the inter- firing time, the equality case 
provides the lower limit of number of targets missed due to 
reloading. For instance, if the reload time unit is twice 
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as long as the inter-flrlng times, a reload requiring five 
time units would cause ten targets to be missed, while if 
the reload time unit is equal to the inter- firing time, only 
five targets would be missed during a reload of five time 
units. As a result, it is assumed that the reload time unit 
is equal to the inter-flrlng time. Consequently, this study 
yields a lower limit on targets missed due to reloading. 

This is the same as saying that this study provides an upper 
limit on engagement capability when faced with a saturation 
attack. 
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IV. THE MODEL 



A. OPERATION PROCEDURES 

Once a raid is detected, the battery commander selects 
a target, decides upon a suitable warhead, and orders the 
launching area to electrically prepare a missile containing 
the desired warhead for firing. At this point, the launching 
area either has the appropriate warhead available or not. 
Whether the desired one of the five possible types is in a 
position to be fired (on a launcher) depends upon the mix 
of the basic load, the order in which the basic load was 
stored, and the reloading policy. As was mentioned previously, 
the mix is assumed to be fixed by some external agency. The 
decision variable available to the battery is the order of 
storage of the warheads (which launcher or what position in 
storage would be the best location for each missile). This 
decision variable together with the reloading policies forms 
the battery's control of its ability to engage targets. 

The question of optimal reloading doctrine is very dif- 
ficult because at least two policies are possible. The first 
is to move only enough missiles to refill empty launchers. 

This allows short reload times at the beginning of the air 
battle. However, as the battle continues, reload times ra- 
pidly increase because the missiles have to be moved farther. 
For instance, suppose two missiles have been fired when a 
section enters reload status. In this case, the first and 



18 



second missiles in storage could be moved onto the launchers 
very quickly. Now suppose storage locations #1 to are 
empty due to previous reloading. When reloading is again 
commenced, the missile in storage #5 (the closest available) 
obviously has to be moved across the four empty storage 
positions. Naturally this takes far longer than moving the 
missile in storage #2 onto a launcher. 

The second possible policy is to move all missiles as 
far as possible each time the section reloads. This policy 
causes the longest reload times to be at the beginning of 
the air battle and hence, causes the reload times to decrease 
with each successive reload. For instance, suppose one mis- 
sile has been fired when reload occurs. The missile in 
storage #1 would be moved onto the launchers, storage #2 
would be moved to storage #1, etc., until all missiles in 
storage have been moved. In total, six missiles would be 
moved one position each. It should be noticed that at the 
finish of reloading, the empty storage location is at the 
far end of the storage barn l.e. storage #6 is empty. No 
missile is ever moved across empty storage locations under 
this policy. 

The contrast between the two policies is easily dis- 
cerned when a simpler system is considered, l.e. only one 
launcher and six storage locations. Under the first policy, 
the first reload takes one time unit (storage #1 moves to the 
launcher), the second reload takes two time units (storage 
#2 moves across storage HI to the launcher) and so on. A 
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total of twenty-one time units would elapse during reloading. 
This is the sum; l+2+3+^+5+6. Under the second policy, 
the first reload takes six time units (storage #1 moves to 
the launcher, storage #2 moves to storage #1, ..., storage 
#6 moves to storage #5). During the second reload, five time 
units would elapse (five missiles move one position) and so 
on.' Thus, total time elapsed is 6+5+^+3+2+l=21 . The dif- 
ference in order of the two equal sums reflects the differ- 
ence in the two policies. Because the first policy allows 
the shortest reload times at the beginning of the battle, it 
is assumed that this reloading policy would be followed. 

This has the intuitive appeal that during any attack of short 
duration the minimum time will be lost due to reload policy 
while during any attack of sufficiently long duration to 
exhaust all missiles in the section a decision maker would 
be indifferent between the two policies. 

It should be noted that this policy assumes that each 
reload was a complete reload. For example, if two launchers 
were empty when the section stopped to reload, both launchers 
were filled. The difference in the two reloading policies 
resulted from movement inside the storage barn. It certainly 
is a feasible policy not to completely reload each time. 
However, the decision rules on when to and when not to reload 
completely are very tenuous and thus are extremely difficult, 
if not impossible, to model. Hence, the model in this study 
includes a complete reloading policy. 
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Additionally, another reloading policy has to be estab- 
lished; this policy Concerns when the section reloads as 
opposed to how the section reloads as discussed above. The 
section has to reload at least twice: fire three missiles, 
reload, fire three missiles, reload. At the most the sec- 
tion would reload six times: fire first missile, reload, 

fire second missile, reload, . . . , fire sixth missile, 
reload. 

A possible policy is to require a reload after a fixed 
number of missiles have been fired, e.g. always reload after 
two missiles have been fired. However, this policy ignores 
the purpose of the section l.e. provide warheads/missiles 
requested by the battery commander to engage the target, 
because the next request might have been for the last mis- 
sile available. Rather than use such a policy, this study 
assumes a policy based upon the series of demands for war- 
heads. The policy is to stop to reload whenever a demand 
cannot be satisfied exactly. This policy Includes the fea- 
ture that if the request for a warhead cannot be satisfied 
exactly, the request is shifted to another section. If the 
request cannot be satisfied in the battery, the target is 
allowed to pass unengaged. This assumption is not as 
reckless as it immediately appears. This type of missile 
unit is considered to be part of a defense-in-depth. Thus 
the fact that one battery does not engage a target does not 
imply that the target reaches its objective unscathed 
because batteries further behind could possibly engage it. 






21 



Of course. If the battery is not a part of a defense-in-' 
depth or is the last battery in line, such a policy would 
be questionable. However, for this study it is assumed that 
this is a realistic procedure. 

The exact response procedure requires some consideration. 
If the warheads can be placed in some order of effectiveness, 
such as size of explosive warhead, it is conceivable that a 
firing policy could exist in which a request for one type of 
warhead could be answered with that exact type or with a 
bigger type warhead. This policy would be a realization of 
a minimum-assurance-of-destruction philosophy. Unfortunately, 
further consideration makes it obvious that under such firing 
policy, the battery would soon be left with only small war- 
heads. Additionally, it is conceivable that the warheads 
cannot be placed in an order of effectiveness. For instance, 
the warheads could be designed for such diverse objectives 
that they could have only a nominal rather than an ordinal 
ordering. As a result of these considerations as well as 
the defense-in-depth concept, it is assumed that the firing 
policy will be an exact response or none. 

Given this exact response procedure, the reason for an 
either fire or reload policy becomes slightly more apparent. 

As was discussed previously, the Interarrival times of tar- 
gets are assumed to be a constant equal to one time unit 
where one time unit is the time required to move a missile 
over one position. Hence, if a request cannot be met, one 
time unit is immediately available for reloading. 
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Consequently, the policy In this study Is to either fire or 
reload depending on whether a desired warhead Is available 
or not. 

Another reloading policy must be established. This 
policy concerns the decision as to which launcher would be 
designated to fill a request when two or more missiles with 
appropriate warheads are on the launchers. For example. If 
launchers #1 and #3 both contain the same type warhead and 
a request arrives for that type warhead, a decision has to 
be made concerning which launcher (#1 or ^3) will be used 
to fill the request. 

Intuitively, It would seem that launcher #3 is the best 
choice because. If the section stops to reload on the next 
request, launcher #3 can be more quickly refilled than If 
#1 Is empty and the missiles on launchers §2 and #3 must 
both be moved over one position and then one missile moved 
out of storage. This procedure Is used In the analysis. 
However, to be sure that Intuition does Indeed hold true. 

Its validity Is checked In the Investigation portion of this 
thesis . 

B. ANALOGY TO A MARKOV CHAIN 

Since the Interarrlval rate of requests Is a constant, 
the stream of requests forms a deterministic queue. As the 
requests arrive, they can either be filled or not filled. 

If It Is filled, a missile Is fired and a counter (the state 
variable) Increases by one. If the request Is not filled, 
reloading occurs and the state variable remains constant 
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for the number of time units required for the reloading 
because firing is impossible during reloading. This con- 
tinues until the state variable reaches a value of nine 
because then the section would be empty. This procedure 
forms a pure birth process with a discrete time space, a 
special case of Markov Chains. The number of time units 
elapsed until the state variable reaches nine indicates the 
number of targets missed due to. reloading because each time 
unit elapsed is equivalent to one target missed since it 
has been assumed that the inter-firing time is equal to the 
reloading time unit . 

A modification of the Markov Chain analogy is necessary. 
Immediately after the finish of reloading, another request 
arrives. If this request cannot be filled, the time variable 
does not increase because no reloading is required (all 
launchers would be filled) and no targets are missed due to 
reloading. As a result, the time parameter increases after 
a reload are conditioned on the arrival of a request that 
could be filled. All intervening unfilled requests are 
ignored. 

Continuing the Markov Chain analogy, this study is 
Interested in the number of time increments elapsing from 
the beginning until the state parameter reaches nine. At 
least nine increments will pass because one missile is fired 
each time increment. As was mentioned above, at least two 
reloads will be necessary and as many as six may occur. 

Hence, the time parameter will be considerably larger than 
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nine. Actually, this study Is concerned with the number of 
time increments during which the state parameter remains 
constant because all possible arrangements of warheads will 
require the nine time Increments to fire the missiles. The 
objective is to discover if it is possible to affect signif- 
icantly the number of time Increments required and its 
resultant of missed targets due to reloading. 



/ 
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V. ANALYSIS 



This analysis will concern a section consisting of 
three launchers and a storage barn for six missiles; thus 
nine missiles will be considered. Additionally, each mis- 
sile is equipped with one of five possible warheads. Thus 
X 2=3 indicates that the missile in position #2 contains a 
type 3 warhead. For the rest of this discussion, the mis- 
sile will be referred to as a type 3 missile since, even 
though the missiles are all constructed identically, the 
different warheads cause each to be used for different 
purposes . 

There exist seven possible configurations of empty 
launchers when the section stops to reload. These seven 
possibilities come from the fact that position #1 or posi- 
tion #2 or position #3 will be empty where the "or" is the 
mathematical inclusive "or". Since there exist three posi- 

o 

tions that are empty or full, there are 2 ~* or eight possible 
configurations of empty launchers. However, one possible 
configuration is that all launchers are full. This possi- 
bility has no physical meaning because the section will not 
stop to reload if it has not fired any missiles. An alter- 
nate and equally valid interpretation is that zero time 
will be required in order to reload. Either way, the con- 
figuration will not occur and hence, 2“^-l = 7 possible 
configurations will occur. 
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Now, the seven possible configurations are not equally 
likely. In fact, the -probability of a certain configuration 
being realized is dependent upon what missiles (that is, 
warheads) are on the launchers and the probability that each 
type of warhead is requested. can be interpreted as the 
probability that a request will be for the warhead on 
launcher #k and equivalently, the probability that launcher 
#k is empty when the section stops to reload. 

For an illustration, the probability that the second 
and third launchers are empty when the section stops to 
reload is 

2y^y (1-y^) 

— =—^ — Cl) 

+ Y2 + Y3 ^ ^ 

The numerator is the probability that two independent events 
occur and that the two can occur in two different orders. 

The factor in parenthesis is the probability of a failure 
on the third request; this is a necessary term to ensure 
that exactly two requests are satisfied. The denominator is 
the conditioning factor. The conditioning is the probability 
that one success has occurred. Hence, the total expression 
is the probability that launchers #2 and #3 are empty when 
the section reloads given that one request has been fired. 
This conditioning is necessary because any targets missed 
before the first firing are Ignored. 

During the reloading process, for this illustration, 
the missile in storage #1 is moved to launcher #2 (2 time 
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units required) and the missile in storage #2 is moved to 
launcher #3(2 time units required). By the model threat 
characteristics described previously, the four time units 
consumed result in four requests unfilled. The section is 
then able to resume answering requests because its launchers 
are filled — but possibly with warheads different from the 
original set; four requests have been missed; and storage 
locations #1 and #2 are empty. The probability that this 
situation is reached when the section is ready for the 
second round is given by Equation (1). The formulae for 
calculating the probability of all seven possible situations 
are given in Table I. 

It must be noted that Equation (1) assumed that all 
warheads on the launchers are different. However, this is 
not certain. The formulae for calculating the probability 
of each, configuration occurring in the four situations in 
which two or three of the three warheads are identical are 
given in Tables II, III, IV, and V. Note that in Tables 
I to V, the probability formulae sum to one. This is 
necessary since each table is an exhaustive list of all 
possible configurations when the section stops to reload. 
Additionally, Tables II to V Indicate that some of the pos- 
sible configurations of empty launchers have a probability 
of zero. This is due to the reloading doctrine requiring 
that when more than one launcher could fill a request, the 
one closest to the storage barns will be used. 
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In order to discover how many requests are not filled, 
it Is only necessary to sum the number of requests that 
are missed each time the section reloads until the storage 
barn is emptied — for then the section would no longer stop 
to reload. The probability of missing that sum would be 
the probability of realizing that particular series of 
empty launcher configurations. 

As an illustration, consider a basic load of warheads 
emplaced in the following manner: {123^5123^} 

where {1 2 3) are on the launchers and {^1512 3^} are in 
the storage barn. Suppose all three launchers are emptied. 
The probability of this event is (from Table I since all 
warheads on launchers are different): 



During this reload, nine time units would be required (stor- 
age missile to launcher #1, ... , storage missile #6 to 
launcher #3) and consequently, nine requests are unfilled. 

The section is then ready to fire again. Suppose all three 
launchers are emptied again. This reloading will cause 
eighteen requests to be unfilled and the probability of this 
happening is Equation (2) again. Hence, the probability of 
twenty-seven (nine plus eighteen) requests being Unfilled is: 




( 2 ) 



yi + Y2 + 




( 3 ) 
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This procedure can be followed for all possible sequences 
of reloading by keeping track of requests missed and multi- 
plying appropriate probabilities as required from Tables 
I to V. 

Unfortunately, there exist a very large number of paths 
in the resulting decision tree. This is evident when it is 
considered that each node can have up to seven paths leading 
from it , It is true that one branch ends after two nodes 
(the one described immediately above), but other branches 
would continue for six nodes (only one missile fired each 
time the section stopped to reload) . Consequently, an 
evaluation of the distribution of targets missed has to be 
accomplished by either computer calculation of the probabil- 
ities in the decision tree or computer simulation of the 

/ 

system. 

In order to have some means of judging the validity of 
the computer calculations (see Computer Program I), three 
criteria are used. The first criterion is that the calcu- 
lated mass function for number of requests not filled must 
sum to one. This is not as trivial as it sounds. The pro- 
gram is designed to feature separate calculations for each 
possible branch and then add the calculated branch mass to 
the mass of the appropriate point in the domain of the missed 
request random variable. If, after all possible branch 
masses have been calculated, the sum is Indeed one, it can 
be felt that the program covers all possible branches. This 
indirect type of assurance of validity is needed for such a 
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program because it Is impossible to manually verify the 
results. Additionally, the computer solution can keep track 
of the number of branches that yield each particular value 
of targets missed. 

Another criterion established to judge the validity of 
the computer solution Involves the minimum possible results. 
Examination of the reloading system reveals a situation 
that appears to yield the minimum number of targets missed. 
The situation is that only one launcher is empty each time 
reloading is commenced and that this launcher is the one 
closest to the storage barn. In other words, six reloads 
are required, each one with launcher #3 empty. The total 
requests missed in this case is twenty-one. In general, 
for a section configured such as the one in this study and 
with this measure of effectiveness, the minimum number of 
failures to supply requested warheads will be: 



M 



min 



S 

E 1 
1=1 



( 4 ) 



The last criterion concerns the maximum number of 
requests missed. This arises when only one missile has been 
fired between each reload but launcher #1 is empty rather 
than launcher #3. In other words. Instead of the empty 
launcher being the closest to storage, the empty space is 
the farthest from storage. The first reload would cause 
three requests to be missed because launcher #2 must be 
moved to launcher #1, launcher #3 then moved to launcher #2, 
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and finally storage #1 moved to launcher § 3 , In similar 
fashion, the second reload causes four requests to be missed 
until the sixth and last reload causes eight requests to be 
missed. 

In summary, the total number of targets missed will be 

8 8 2 

Zi = Zi-Zi 
i=3 i=l i=l 

The right hand side has the intuitive content that the 
first term predicts the mamber of unfilled requests if the 
system consisted of one launcher and a storage of eight. 

This is completely analogous to Equation (4). But since 
this was not a single laimcher system, the second term is 
a correction term for the intervening launchers . The second 
term is the number of requests missed if the missiles on 
launchers #2 and #3 would indeed have to be moved to 
launcher #1. Since these two missiles do not move in such 
fashion, the sum is subtracted from the total. This result 
is easily expanded to the general case: 



M 



max 



S+L-1 
Z 1 
1=1 



L-1 
Z 1 
1=1 



( 6 ) 



The analytical program was run on several test cases 
utilizing Tables I to V. In all cases the calculated mass 
function summed to a number larger than .9999 and had a 
domain from 21 to 33 (the limits predicted by Equations (4) 
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and (6)). Hence, it appeared that the program was indeed 
performing properly and thus could be used to calculate the 
mass functions of number of targets for any possible 
permutation of warheads and probabilities of requests for 
each type of warhead. 

During conduct of this trial phase, one surprising 
result was noticed. The minimum number of targets missed 
in all cases was twenty-one as predicted. However, the 
output of the program indicated that this result was 
achieved in as many as four ways Instead of only one way 
(the one used to establish Equation (4)). It was not pos- 
sible to easily find the alternate paths to this minimum 
result because as many as 4529 paths existed in some cases. 
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VI. SIMULATION 



A method of validating computer solutions is to find 
the solution by another method. In this situation, computer 
simulation is a convenient alternative. Accordingly, a 
simulation was constructed embodying all the characteristics 
discussed previously (see Computer Program II). A random 
number generator was Included to produce a series of 
requests of warheads in accordance with the Pj . With the 
request, the computer then determined if a matching warhead 
was on the launchers. If so, the position (Xj^) was set to 
a value of zero and another request was generated. If not, 
reloading was accomplished, a record of requests missed 
during reloading was updated, and then the next request was 
generated. This procedure continued until the storage barn 
was emptied. Then a counter corresponding to the number of 
requests missed was Increased (if 27 requests were missed, 
the counter for event "27" was Increased by one), the 
section was reinitialized to the original basic load of 
warheads, and a new stream of requests was presented to the 
section. This procedure was performed for a fixed number 
of repetitions and then an observed mass function was cal- 
culated by dividing the observed occurrence of each event, 
e.g. total number of times that 23 requests v;ere missed, by 
the total number of repetitions. Additionally, records were 
maintained of the requests in order to ensure the proper 
proportion of types of requests were Indeed being presented 
to the section. 
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Table VI presents the simulated mass function of the 
random variable, requests missed, for varying numbers of 
repetitions. It will be noticed that a very large number 
of repetitions was necessary In order to produce a steady 
mass function. 

Using the K.Olmogorov-Smlrnoff test to determine If the 
simulation results from a sample of 15,000 could have been 
realized from the analytical mass function [Siegel 1956], 

It was determined that the mass function yielded by the 
simulation did Indeed come from the analytical mass function 
that Is, the two were Identically distributed. The size of 
this hypothesis test was .05 [Zehna 1970]. In other words, 
there existed a 5^ probability that the test could have 
declared the two mass functions different when In fact they 
were Identical (See Table VII). This result was very encour 
aging since the author had very high confidence that the 
simulation accurately performed In the manner of the modeled 
system. With the hypothesis test revealing that the simula- 
tion and the analytical model produced the same results, it 
was felt that the analytical work was accurate. 
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VII. INVESTIGATION OP THE SYSTEM'S CHARACTERISTICS 



Since the number of requests missed due to reloading 
time is a random variable, some means are necessary to 
determine if one mass function is preferable to another 
mass function. The different mass functions would result 
from different ordering of the warheads in the section. 

The most familiar criterion and easiest to calculate and 
understand for ranking desirability of random variables is 
that of relative sizes of expected values. This criterion 
ranks mass functions in order of their expected values. In 
this study, the objective is to minimize requests missed due 
to reloading. Hence, the arrangement of warheads that 
produces the smallest expected value of the random variable 
would be the most desirable arrangement. 

However, another more inclusive criterion is available. 
This is ordering by survivor functions. For a random vari- 
able M, the survivor function is defined as: 

F(m) E Prob{M>m} ( 7 ) 

where m is an element of the domain of the random variable M. 
Note that the survivor function satisfies 

F(m) = 1.0 - P(m) (8) 

where F(m) is the well known cumulative distribution function. 
For two random variables and M2 (resulting from two war- 
head arrangements), it is true that if F^(t) ^ F2(t) then 
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E(lVIj^) ^ E(M 2 ) . Hence, If while comparing two particular 
warhead arrangements. It could be shown that P-j^(t) ^ P 2 (t), 
a considerably stronger ranking would be known because the 
expected value ranking would be Identical. On the other 
hand, if it is shown that the two expected values are dif- 
ferent, nothing is known (without further calculations) 
about the stronger ranking. 

With this knowledge, it can then be determined what 
warhead orderings are preferable to others. Por example. 
Tables VIII, IX, and X provide the expected values and sur- 
vivor functions for three different arrangements of one mix 
of warheads. Examination reveals that E^g < E^ < Eg with E^ 
and Eg virtually equal. Also it is noted that F^g £ P^ and 
P^g '£ Pg while Pg and Pg have no clear cut dominance. This 
can be expected since their expected values are so close to 
each other. Thus the warhead arrangement in Table X is 
preferable because its survivor function dominates the other 
two. 

Examination of P^ reveals a discrepancy. P^ does not 
become .0000 at m=3^ as necessary because m=33 was the 
largest value of m with a nonzero mass. The explanation is 
that the analytic program yields a mass function that does 
not sum exactly to one and the difference appears as the 
nonzero value of P^ at m=3^. 

At this point, it is apparent how the best warhead 
arrangement can be determined l.e. find the permutation of 
warheads that yields a mass function that dominates all 
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other possible mass functions. Unfortunately, there exist 
numerous possible permutations of any given mix of nine 
warheads. The upper limit of possible permutations would 
occur when the mix is 2:2: 2: 2:1. In this case, the number 
of permutations would be 

9 ! 

( 2 !)( 2 !)( 2 !)( 2 !)( 1 !) " 

[Hickman 1971]. It is Immediately evident that the size of 
this task is overwhelming. If it was judged worthwhile, the 
optimal permutation could be obtained by using the analytical 
computer program to calculate the mass function for every one 
of the possible permutations and then use the two criteria 
described above to find the best permutation. It is cer- 
tainly conceivable that a unique solution would not exist. 
However, the criteria would still be able to eliminate poor 
permutations until a set was delineated that could not be 
differentiated further. At this point, an arbitrary selec- 
tion would be in order. It is noted that this process in 
essence eliminates poor choices and then selects a solution 
from the remaining possibilities instead of finding the 
optimal solution directly. 

Since finding the set of "best" permutations involves 
evaluating all possible permutations, some faster method 
is desirable because the computer analytical evaluation of 
each mass function requires approximately four seconds. 

The upper limit on possible run times is then slightly under 
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eight and one-half hours! In an effort to determine char- 
acteristics of the arrangement in Table X which were not 
present in the arrangements in Tables VIII and IX, it was 
noticed that the warheads were placed in a nondecreasing 
order of likelihood of request, e.g. p^ warheads were on 
the launchers while p^ warheads were in the rear of the 

storage barn when P2^^?2^P3^P4^P5 • 

As a result, a factorial computer experiment was con- 
ducted. Two effects were considered: (1) the effects of 

the monotone arrangement of warheads, and (2) the effects 
of the firing policy of firing the launcher closest to the 
barn when two missiles on the launchers contained the 
desired warhead. The question concerning the latter effect 
was explained earlier. The experiment was conducted with 
the simulation since it had been shown to yield very accu- 
rate results and was easier to modify than the analytic 
program. 

Comparison of parts A and B of Table XI revealed that 
the firing doctrine had virtually no effect on the perform- 
ance of the section loaded in the manner shown. This seemed 
reasonable since the doctrine concerned situations in which 
two or three identical warheads were on the launchers. With 
the loading plan as shown, duplication of warheads would not 
occur very often and hence, the firing doctrine would have 
little opportunity to manifest Itself. 

Examination of Tables XII and XIII revealed that, for 
the permutations of warheads which caused frequent duplica- 
tion of warheads on the launchers, the firing doctrine 
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under study resulted in a dominance by expected values and 
by survivor functions in both tables. Hence, the assumption 
previously made to use this doctrine appears to be 
reasonable . 

With this doctrine established as beneficial, examination 
of part A of Tables XI, XII, and XIII revealed that Table 
XIII dominated both Tables XI and XII at the survivor func- 
tion level. It was also indicated that no particular 
dominance existed between Tables XI and XII. Even though 
the expected value of Table XI was less than that in Table 
XII, the difference was so slight that any declaration of 
dominance would be reckless because simulation was used to 
obtain the results. 

Af'this point it was again noticeable that the permuta- 
tion of warheads with a nondecreasing ordering of warheads 
yielded the dominant mass function. Many more trials could 
have been conducted, but since so many orderings were 
possible (perhaps, it was not necessary for a nondecreasing 
ordering but merely an overall upward trend in probability) 
the procedure seemed to promise little additional knowledge 
in return for much more computer run time. In fact, several 
more trials were conducted with no results contradicting the 
principle delineated above. 

Instead, an explanation was discovered for why the non- 
decreasing orderings appeared to be the rule of dominant 
permutations. Based upon the probability formulae in Tables 
I to V and the results of considering the firing doctrine 



of closest-launcher versus farthest-launcher in case of 
duplication, it was discernable that, in terms of requests 
missed during reload, it was preferable to have the empty 
launchers closest to the storage barn whenever possible. 
Actually this is in keeping with judgement: if two launchers 

are empty, it is best that they are launchers #2 and #3 
because missiles can be moved directly from storage to a 
launcher without losing time moving a missile between 
launchers. Subsequently, the question arose: given three 

warheads to be placed on the launchers, what is the best 
arrangement to ensure that launcher #3 would be empty more 
often than #2 and that #2 would be empty more often than #1? 
If the three warheads are Identical, the firing doctrine 
could ensure the desired configuration; if the three war- 
heads are not identical, it is obvious that the best arrange- 
ment is to put the warhead most likely to be requested on 
launcher #3» the warhead next most likely to be requested on 
launcher #2, etc. 

When consideration is given to what warhead should be 

placed in storage #1 after the three on the launchers have 

been placed, it appears that the same principle applies. 

The warhead in storage #1 should have a p. greater than or 

J 

equal to y^ because the warheads on the launchers have 
already been arranged so that launcher #3 would most likely 
be empty. With launcher #3 empty, the missile in storage #1 
would be moved to the launcher. When this occurs, it would 
still be desired that y^ 3. ^2 — ^1' 



This would be 



accomplished when storage #1 has a pj greater than or equal 
to y^. This reasoning is then immediately extended to all 
positions in the storage barn;. In conclusion, it can be 
established that the number of requests not filled due to 
reloading will be minimized by arranging the given assortment 
of warheads in an order beginning with the least requested 
warhead on launcher #1, an equally or more likely warhead on 
launcher #2, and so on until all' nine warheads are positioned. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



At this point, the object of the study has been 
achieved: Specifically, it is known what principle of 

ordering warheads in a sequential storage system will 
result in the minimum number of unfilled requests. However, 
an examination of the principle reveals a startling conclu- 
sion: the best ordering for minimizing misses due to 

reloading causes the worst performance in terms of availa- 
bility! This is immediately evident when the permutation 
in Table X is considered. With the reloading policies as 
specified in this study, no penalty is suffered for requests 
missed until the first success. The permutation enhanced 
the probability that this first success would be on launcher 
#3» a situation shown to be best. Then since the least 
likely requested warheads were on the launchers, it is most 
likely that the next request will not be filled and by 
policy, the section then stops to reload. It was indicated 
earlier that the sequence of firing once and then reloading 
led to the minimum number of misses during reloading, i.e. 

21 misses. 

But consideration has to be given to the number of 
misses until the first success after each reloading. This 
number is a random variable with a geometric distribution. 
The single parameter of this distribution is the probability 
of a success on a single request. In this case, the 
parameter is the probability that a request is for a warhead 
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on the launchers. However, the optimal ordering discovered 
earlier placed the least probable warheads on the launchers. 
Hence, the size of the parameter would be small. The con- 
sequence Is that, since the expected value of a geometric 
distribution Is the reciprocal of the parameter, the number 
of misses before the first success would be large. In 
summary, the ordering of warheads that minimizes number of 
requests missed due to reloading simultaneously maximizes 
the number of requests missed due to availability! The 
Immediate conclusion Is that the optimal ordering sequence 
Is extremely sensitive to the objective function. 

More formally, consider an objective function: 

T 

M = E (a. + m. ) + a^ . (10) 

s 1^1 1 1 

This form of the objective function Is more meaningful when 
the sequence of events In a section Is considered. When the 
section begins the battle, a number of requests will be 
missed at first (a^) because of the type of warheads available 
on the launchers. This number would be a geometric random 
variable whose parameter depends on y^, 3^3* Then 

either one, two, or three requests will be filled after which 
the section will stop to reload. A certain number of 
requests would be missed during reloading (m^). This number 
would be a random variable depending on which launchers were 
empty and which storage locations were filled. After 
reloading, requests would again be missed due to availability 



{ 0 ^ 2 ^ ' This variable would again be distributed geometrically 
but with a parameter different from the first. At least one 
success will be experienced and then reloading will again 
occur with m 2 requests missed. This will continue until the 
last reload (T) empties the storage barn. T is a discrete 
random variable with a domain of 2, 3, 5> 6 (as was 

explained earlier) . The mass function of this discrete 
random variable is not readily apparent. Finally, after the 
storage barn is empty, a certain number of requests will be 

f 

missed until all warheads on the launchers are requested (a ) . 
This number will be similar to a negative binomial random 
variable with a parameter depending on the p. of the last 

J 

three warheads . 

Rearrangement of (10) yields 

T T . 

= Z a. + E m, + a^ . (11) 

® 1=1 ^ 1=1 ^ 

where the first term is a random sum of a variety of possible 
geometric distributions, and the second term is a random sum 
of random misses and the last term is similar but not iden- 
tical to a negative binomial distribution. Hence, M is the 
sum of three random variables each with a rather complex 
mass function. As a result, the mass function of M would be 
extremely difficult to develop analytically. 

This study concentrated on the second term by setting 
the first and third terms to zero by the defense-ln-depth 
capability. The analysis shown on Tables I to V provided 



^5 



the mass function of the second term. The result was that 
the objective of this study was to minimize 

T 

B_ = E m. . (12) 

This was found to be accomplished by placing the' least 
likely warheads on the launchers first. Examination of 
this result revealed that this policy would make 

f 

A = Z a. + a^ (13) 

® i=l ^ 

a comparatively large number. Intuition based upon earlier 
comments indicated that if the objective had been to minimize 
A , the optimal policy would be the opposite of that found 

o 

in the previous chapter, i.e. most likely warheads placed 
on the launchers and each warhead’s location based on 
descending size of Pj . Further intuition also indicates 
that if the objective had been to minimize M^, that is mini- 
mize the number of targets missed due to both reloading and 
availability, the optimal arrangement would be a compromise 
between the two extreme arrangements called for by minimizing 
either A or B_ individually. Of course, these two possible 
objectives should be subjected to formal analysis in order 
to determine if the insights gained in this thesis were 
sufficient for intuition to Indeed provide the correct 
solutions . 
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After that is accomplished, investigation could be 
conducted on the optimal ordering for a battery. In this 
case 




m. 



( 1^0 



since a battery consists of three sections. Substitution of 
(11) into (l4) yields 



®B 



3 T 3 f 

I I a^j + Z a. 

J=1 1=1 j=l J j 



3 T 

+ Z Z m . j 
j = l 1=1 



( 15 ) 



Analogous to the section, it could be assumed that 
missing a request due to availability does not incur any 
penalty because other batteries exist to handle targets 
missed. In which case, the objective would be to minimize 
Bg. Similarly, the objective could be to minimize Ag or to 
minimize Mg. At this level, suboptimization might also 
arise in that the policy that optimizes one section might 
not optimize the battery if all three sections were loaded 
in accordance with an optimal section plan. For example, as 
was mentioned above, a section ordering plan to minimize M 

o 

would likely be a compromise between the extreme strategies 

demanded to minimize A or B . However, the optimal ordering 

s s 

for a Mg might consist of two sections ordered to minimize A^ 
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and one section ordered to minimize B„ rather than all 

o 

three ordered to minimize . 

s 

In a continuation of the structuring effort, an objective 
for an entire defense would be 




B 

Z M, 
K=1 ^ 
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Z 



k=l 




T 

I 

1=1 




B 3 
z • z 

k=l j=l 




( 16 ) 



B 3 T 

+ ^ ^ ^ mkil 

k=l j=l i=l 



At this level of aggregation all terms must be considered 
for there are no batteries behind the defense. Thus targets 
missed due to both availability and reload time must be 
considered. To reflect the relative importance of targets 
missed as the position of a battery in a defense varies, a 
weighting factor might be Introduced. For instance, if a 
defense of five batteries is organized with two batteries 
on line followed by a third battery followed by the fourth 
and fifth batteries on line, targets missed due to availa- 
bility by batteries #1 and #2 would have little consequence 
because a back up capability exists. However, batteries 
and §5 must stop all targets. Hence, let qj^ denote a 
weighting factor of a target missed by battery k. For 
example, a possible set of weighting factors is: 
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1 



Ql = Q2 = 

Q3 = 3 (17) 

Qlj = = 10 



These reflect rapidly increasing importance to a target 
being missed as it penetrates the defense. 

Using this weighting factor in Equation (16), the 
objective function to be minimized is: 



M 



D 



B 3 T 
E E E 
k=l j=l i=l 






+ 



B 

E 

k=l 



3 

E 

j=l 



^k^kj 



(18) 



where q^^ and Sj^ allow for possibly different weighting 
factors for targets missed due to availability and for 
targets missed due to reloading. 

Again it appeared that suboptimization could present a 
problem at this defense level. The optimal defense ordering 
of warheads may not be five (or however many batteries there 
are) optimal battery arrangements. In fact, if weighting 
factors are used, the optimal defense ordering policy almost 
certainly would not be five optimal battery arrangements. 

Another possible form of an objective function Involves 
a slightly different weighting value than those in (17)* 

This method further places a weight upon what type of 
targets are missed by which battery. It was assumed that 
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five types of targets were possible. Hence, the objective 
function to be minimized would be: 




B 3 T 

k=i j=i 1^1 ^^kAji ®ia\ji^ 



B 



+ ■2 2 
k=l jfl ^ 



(19) 



In this manner, the penalty for battery k missing a type ^ 
target due to availability may be different from the same 
battery missing the same target due to reloading. Using 
misses due to reloading as an example, the weighting factor 
for a battery in the first echelon of defense missing a 
nuclear weapon carrier would be less than the weighting 
factor for a battery in the last echelon missing the same 
target. On the other hand, one of the five possible types 
of targets may be of such low priority that no difference 
in weighting factors between batteries is warranted. This 
last objective function provides enough flexibility for 
weightings for different missions and goals that virtually 
any situation could be modeled. Unfortunately, solutions 
of these possible objective functions would be quite 
complicated. However, based upon the results of the section 
analysis. It Is apparent that suboptimization is a very real 
possibility in this type of system and consequently, any 
further analysis must be conducted at the defense-level if 
possible . 
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TABLE I: 



Probability Formulae for 

y^¥y^ 12 2 3 



Position of 

Empty Launcher Probability of Occurrence 

12 3 



X 



X 



X 



X X 



X X 



X X 



XXX 



YlCl - Y 2 

y]_ + y2 + y3 



y2(i - yi - ^3) 

yi + Yp + y3 



y^Ci - y^ - y2) 

yi + y2 + Y3 



Py^y^Cl - y^) 
y-^ + y~2 + y3 



2y^y2(l -y^) 

Yl +Y 2 +Y 3 



2y2Y3(i - y^) 



Yl + Y 2 + Y 3 

6y3^y2Y3 

Yl + Y2 + y3 
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TABLE II: Probability Formulae for y2~y3» y^^y2 



Position of 

Empty Launcher Probability of Occurrence 

12 3 



y2(l - y;^) 






- yj) 



zyiygd-y^) 

yi+yj 



3 y 2 y 2 

y^TyJ 
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TABLE III: 


Probability Formulae for y2=y^ y^^?^y 2 


Position of 
Empty Launcher 


Probability of Occurrence 


12 3 




X 


y^^d - ^2^ 

yi + yj ■ 


X 


0 


X 


y2d " yi " y2^ 

yi + y2 


X X 


0 


X X 


2y^y2(l - ^2^ 
yi + yj 


X X 


yj^d - y^) 
y^ + ys 


XXX 


3yj^y2^ 

yi + yj 
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TABLE IV: Probability Formulae for 



Position of 
Empty Launcher 


Probability of Occurrence 


12 3 




X 


0 


X 


Yid - - y3) 

yi + yj 


X 


y^(l- y^) 
yi+ys 


X X 


yj^^Cl - yj) 
yi + ys 


X X 


0 


X X 


2y]_y3(i - y]_) 

yi + yj 


XXX 


3yj^^y3 
yi + yj 
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TABLE V; Probability Formulae for y^^=y2=y2 



Position of 

Empty Launcher Probability of Occurrence 

12 3 



X X 



X X 



XXX 



1 - yi 



yid - y^^) 



55 



TABLE VI: Results of Simulation with Various Number of Replications 



o 

o 

o 

LTV 

I — I 



OO^VOOOrHOOOOOC^mOJLTNiHO 
o o o ^ CO ON oo^ c\j m o o o co 

O O O O MD iH OO LTN CTN cr\C\J O O O O OO 

oooooiHojojiHoooooo cn 



m 



LO 

o 



o 

o 

o 



o o ^ oj v£) m OO cTNJ^r o o oooo m o 

OOOirNLTSOMD OLr^^-OOLr^OOO 
OOOOMDOJOOLncr\C^C\JOOOO 
OOOOOiHCUOJiHOOOOOO 



o 

VO 



OO 

C\J 



o 

OJ 



m 

OJ 



o 

o 

o 

LTN 



OOOO^OOOOJ^^OJ^OOO 
O OOLnLnOJLTNiHVOmOOLrNOJOO 

oooovdojooltnctnc^ojoooo 

OOOOOiHOJOJiHOOOOOO 



<x> 



1-3 1-3 



m 



iT\ 



^ =T 



o 

o 

o 



000000000000000 

OOOOLTNOOOOt^iHCTN^iHOO iH 
OOOiHVD^OJLTNOOOiHOOOO OJ 
OOOOOiHOJOJiHiHOOOOO 



CTn 

00 



00 

C\J 



VO 



o 

o 

LTN 



000000000000000 
ooooooojvoojoooo-=rooo 00 

O O O iH LTN VO CTnCTnOJ O O O O iH 

OOOOOiHOJOJiHOOOOOO 



00 

OJ 



00 



OJ 



o 

o 



000000000000000 
000000000000000 
OOOH^OD— ^^t>-moooo 
OOOOOHOJOJOJOOOOOO 




o H c\j ooj=r mvD t>-oo ctno h oj 00^ 
cvjojojojojojojojojojoooooooocn 



I 



g 

S' 



0 ) 
O CO 
0 ) 

^5 
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(each column Is a simulated mass function) 



TABLE VII: Kolmogorov-Smlrnof f Test Validating the 

Simulation in One Case 



H : the mass function of the simulation is the same as the 

analytical mass function 

H,: the mass function of the simulation Is different from 

the analytical mass function 



o = .05 



crlt 



1.36 



Vl5,000 



= .0111 



12 3 ^ 5 6 7 8 9 
5 i 2 3 ^ 



1 2 3 



3^5 



Pj .ilO .25 .20 .10 .05 



X 


Fo(x) 


^15,000^^^ 


-^15, 


20 


.00000 


.00000 • 


.00000 


21 


.00000 


.00000 


.00000 


22 


.00032 


.00040 


.00008 


23 


.00^17 


.00500 


.00083 


2H 


.072^1 


.07327 


.00086 


25 


.19041 


.19240 


.00199 


26 


.42086 


.42573 


.00487 


27 


.66976 


.67973 


.00997 


28 


.86760 


.87206 


.00446 


29 


.96676 


.96992 


.00316 


30 


.99409 


.99419 


.00010 


31 


.99925 


.99939 


.00014 


32 


.99987 


.99992 


.00015 


33 


.99990 


.99999 


.00009 


3^ 


1.00000 


1.00000 


.00000 


Since d„_^ 
max 


= .00997 < 


.0111 = d , accept . 

crlt 0 





(x) 
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TABLE VIII: Analytical Mass and Survivor Functions for 

a Permutation of Warheads with Decreasing 
Order of Probability of Request 



12345 6 7 8 9 



1 1 



3 3 ^^ 4 5 



J 12345 
.40 .25 .20 .10 .05 



Number of 
Targets Missed 



Mass 



Survivor 
Function (F0) 



20 


.0000 


1.0000 


21 


.0001 


.9999 


22 


.0016 


. .9983 


23 


.0049 


.9935 


24 


.0403 


.9532 


25 


.1573 


.7959 


26 


.1772 


.6187 


27 


.3062 


.3125 


28 


.1747 


.1378 


29 


.1097 


.0281 


30 


.0281 


.0000 


31 


.0000 


.0000 


32 


.0000 


.0000 


33 


.0000 


.0000 


34 


.0000 





Expected value (E0) = 26.837 
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TABLE IX: 



Analytical Mass and Survivor Functions for a 
Permutation of Warheads with an Arbitrary 
Order of Probability of Request 



1 2 3 ^ 5 6789 



3 ^ 5 



3 ^ 



j 123^5 
p. .ilO .25 .20 .10 .05 



Number of 
Targets Missed 



Mass 



Survivor _ 
Function (F^) 



20 


.0000 


1.0000 


21 


.0000 


1.0000 


22 


.0003 


• .9997 


23 


.0038 


.9958 




.0682 


.9276 


25 


.1180 


.8096 


26 


.2304 


.5791 


27 


.2489 


.3302 


28 


.1978 


.1324 


29 


.0992 


.0332 


30 


.0273 


.0059 


31 


.0052 


.0007 


32 


.0006 


.0001 


33 


.0000 


.0001 


3 ^ 


.0000 


.0001 



Expected value (E^) = 26.811 
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TABLE X: Analytical Mass and Survivor Functions for a 

Permutation of Warheads with an Increasing 
Order of Probability of Request 


1 12 3 


456 


789 


J 12345 


XjL 5 4 4 


332 


2 11 


Pj .40 .25 .20 .10 .05 


Number of 
Targets Missed 




Mass 


Survivor 
Function (^]^q) 


20 




.0000 


1.0000 


• 21 




.0318 


.9682 


22 




.0843 


.8838 


23 




.2156 


.6682 


24 




.2125 


.4557 


25 




.1818 


.2739 


.26 




‘.1433 


.1306 


27 




.0868 


.0438 


28 




.0326 


.0112 . 


29 




.0104 


.0008 


30 




.0008 


.0000 


31 




.0000 


.0000 


32 




.0000 


.0000 


33 




.0000 


.0000 


34 




.0000 


.0000 



Expected value (E-j^q) = 24.435 
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TABLE XI 



: Comparison of Close-Launcher versus Far- 

Launcher Firing Doctrine for a Permutation 
of Warheads with an Arbitrary Order of 
Probability of Request 



1 123^56789 J 12345 

1 2 3 4 5 i 2 3 4“ Pj .40 .25 .20 .10 .05 



Part A; Close Launcher 



Number of 
Targets Missed 


Mass 


Survivor 

Function 


20 


.0000 


1.0000 


21 


.0000 


1.0000 


22 


.0004 


.9996 


.23 


.0052 


.9944 


24 


.0656 


.9288 


25 


.1205 


.8083 


26 


.2363 


.5720 


27 


.2509 


.3211 


28 


.1954 


.1257 


29 


. 0970 


.0287 


30 


.0230 


.0057 


31 


.0053 


.0004 


32 


.0003 


.0001 


33 


.0001 


.0000 


34 


.0000 


.0000 



Expected value = 26.785 
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TABLE XI (continued) 



Part B: Far-Launcher 



Number of 
Targets Missed 


Mass 


Survivor 

Function 


20 


.0000 


1.0000 


21 


.0000 


1.0000 


22 


.0000 


1.0000 


23 


.0002 


.9998 


24 


.0501 


.9497 


25 


.1026 


.8471 


26 


.2406 


.6065 


27 


.2767 


.3298 


.28 


.2033 


.1265 


29 


.0978 


■.0287 


30 


.0230 


.0057 


31 


.0053 


.0004 


32 


.0003 


.0001 


33 


.0001 


.0000 


34 


.0000 


.0000 



Expected value = 26.89^ 
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TABLE XII; 



Comparison of Close-Launcher versus 
Far-Launcher Firing Doctrine for a 
Permutation of Warheads with a 
Decreasing Order of Probability of 
Request 



i 123456789 J 12345 

I I 2 2 3 3 4 4 Pj .40 .25 .20 .10 .05 



Part A: Close-Launcher 



Number of 
Targets Missed 


Mass 


Survivor 

Function 


20 


.0000 


1.0000 


‘21 


.0001 


.9999 


22 


.0017 


.9982 


23 


.0068 


.9914 


24 


, .0418 


.9496 


25 


.1565 


.7931 


26 


.1712 


.6219 


27 


.3156 


.3063 


28 


.1699 


.1364 


29 


.1107 


.0257 


30 


.0257 


.0000 


31 


.0000 


.0000 


32 


.0000 


.0000 


33 


. ,0000 


.0000 


34 


.0000 


.0000 



Expected value = 26.822 
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TABLE XII (continued) 



Part B: Par-Launcher 



Number of 
Targets Missed 


Mass 


Survivor 

Function 


20 


.0000 


1.0000 


21 


.0001 


.9999 


22 


.0004 


.9995 


23 


.0005 


.9990 


24 


.0149 


.9841 


25 


.0262 


.9579 


26 


.1014 


.8565 


27 


.1599 


.6966 


28 


.2489 


.4477 


‘29 


.2069 


.2408 


30 


.1509 


.0899 


31 


.0642 


.0257 


32 


.0234 


.0023 


33 


.0023 


.0000 


34 


.0000 


.0000 



Expected value = 28.300 
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TABLE XIII: 



Comparison of Close-Launcher versus 
Far-Launcher Firing Doctrine for a 
Permutation of Warheads with an Increasing 
Order of Probability Request 



1123456789 J 12345 
xj 5 4 4 3 3 2 2 I ~ pj . 40 . 25 . 20 . 10.05 



Part A: Close-Launcher 



Number of 
Targets Missed 


Mass 


Survivor 

Function 


20 


.0000 


1.0000 


21 


.0309 


.9691 


22 


.0848 


.8843 


•23 


.2146 


.6697 


24 


.2093 


.4604 


25 


.1803 


.2801 


26 


.1467 


.1334 


27 


.0896 


.0438 


28 


.0318 


.0120 


29 


.0113 


.0007 


30 


.0007 


, .0000 


31 


.0000 


.0000 


32 


.0000 


.0000 


33 


.0000 


.0000 


34 


.0000 


.0000 



Expected value = 24.435 
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TABLE XIII: (Continued) 



Part B: Par-Launcher 



Number of 
Targets Missed 


Mass 


Survivor 

Function 


20 


.0000 


1.0000 


21 


.0000 


1.0000 


22 


.0309 


.9691 


23 


.0932 


.8759 


24 


.1356 


.7403 


25 


.1835 


.5568 


26 


.1974 


.3594 


27 


.1503 


.2091 


•28 


.1128 


.0963 


29 


.0610 


.0353 


30 


.0239 


.0114 


31 


.0100 


.0014 


32 


.0014 


.0000 


33 


.0000 


.0000 


34 


.0000 


.0000 



Expected value = 25.855 



oooo oooo ooo oooo oooo oooo oooo ooooo ooo 



COMPUTER PROGRAM I 



DIMENSION THREAT (5) ,MISSIL( 9), FACTO) ,STOR A ( 27C0 , 1 1 ) , S 
CTCP^B(2700,11) ,PR0B (40 ) , MOVE S (40 ) 

DATA STORA/29700«0.0/,STQRB/29700*0.0/,PRCB/40*0.0/,MO 
CVESMO’f'O/ 

STGRB(1,10)=1.0 

READ IN THE INITIAL ARRAY OF NINE WARHEADS. 

READ (5,8010) (STORBd, J) ,J = 1,9) 

8010 FORMAT (9F2.0) 

READ IN THE PROBABILITIES OF THE FIVE TYPES CF 
WARHEAD REQUESTS. THE FIVE PROBABILITIES MUST ADD TO 
ONE. 

READ (5,8030) THREAT 
8030 FORMAT (5F4.3) 

THIS OUTPUT ENSURES THAT THE INITIAL WARHEAD ARRAY 
WAS READ IN PROPERLY. 

WRITE (6,3040) ( STORB ( 1 , J ) , J= 1 , 9 ) 

8040 FORMAT (‘O', ‘FIRST ARRAY IS‘,9F3.0) 

THIS OUTPUT ENSURES THAT THE WARHEAD REQUEST 
PROBABILITIES WERE INPUT CORRECTLY. 

WRITE (6,8050) THREAT 

8050 FCRMAT ( • 0 • , • WARH EAD PROBABILITIES ARE',5F1C.3) 

THIS STEP DETERMINES WHEN ALL BRANCHES OF THE 
DECISION TREE HAVE BEEN DELINEATED. 

2000 IF (STORBd, 9) .EQ. 0.0) GO TO 6000 

THIS STEP LOADS THE RESULTING ARRAYS FROM THE 
PREVIOUS CALCULATIONS INTO THE WORKING AREA. 

DC 3 JJK=1,2700 
DO 2 JJJ=1,11 

STCRA( JJK, JJJ)=STGRB( JJK, JJJ) 

2 CONTINUE 

3 CONTINUE 

THIS STEP RESETS THE STORAGE AREA TO ZEROS. 

DC 5 JJI=1,2700 
DO 4 JJL=1,11 
STGR6( JJI,JJL)=0.0 

4 CONTINUE 

5 CONTINUE 
NE = 1 

DC 5000 NNN=1,2700 

THIS STEP DETERMINES IF ALL ARRAYS IN THE PRESENT 
CALCULATIONS HAVE BEEN PROCESSED. 

IF (ST0RA(NNN,1) .EQ.0.0) GO TO 2000 

TH-IS STEP DETERMINES THE PROBABILITY FACTORS 
NECESSARY FOR THE CALCULATIONS DICTATED BY TABLES 
I THROUGH V. 

DC 30 J=l,.3 
DC 20 K=l,5 
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noon 



IF (STORA(NNN,J) .EQ.FLOAT(K) ) F ACT ( J ) =THR EAT( K) 

20 CCNTINUE 
3C CONTINUE 
DC AO 1=1,9 
MISSIL(I)=STORA(NNN,I ) 

AC CONTINUE 

IF UMISSIL(1).EQ.MISSIL(2) ) .AND . (MISS I L ( 2 ) . EC.MI SSI L( 
C3)l) GO TO 500 

IF IMISSIL(l) .EQ.MISSIL(2) ) GO TO 390 
IF CMISSIL(2) .EQ.MISSILO) ) GO TO 270 
IF IMISSILd) .EQ.MISSILO) ) GO TO 150 

THE FORMULAE IN TABLE I ARE ENCODED BETWEEN HERE 
AND LINE 150. 

DEN=FACT(1)+FACT(2)+FACT(3) 

MISSIL( 1)=0 
NCRIT=0 

CALL RELOAD (MISSIL,NCRIT) 

IF (MISSIL(9).EQ.O) GO TO 53 
DC 51 N=l,9 
STORB(NB,N)=MISSIL(N) 

51 CCNTINUE 

STGRB(NB,10) = FACT( 1)=«'(1.0-FACT(3}-FACT( 2) ) *ST0RA( NNN , 1 
C0)/DEN 

STCRB(Ne,ll)=STORA(NNN, 11 )+NCRIT 

NE=NB+1 

GO TO 55 

53 ILCS=STGRA(NNN,11)+NCRIT 

SFF.OB=FACT(l)*(l .0-FACT(3)-FACT(2) )*ST0RA(NNN,10) /OEN 
MCVES( ILDS )=M0VES( ILDS)+1 
PRCB(ILCS)=PRO8(ILCS)+SPR0B 
55 DC 57 K=l,9 

M1SSIL(K)=ST0RA(NNN,K) 

57 CCNTINUE 
MISSIL( 2)=0 
NCRIT=0 

CALL RELCAO ( M ISS I L , NC R IT ) 

IF (MISSIL(9) .EQ.O ) GO TO 61 
00 59 M=l,9 
STGR3(NB,MJ=MISSIL(M) 

59 CCNTINUE 

STORB(NB,10)=FACT (2)*(1.0-FACT( 1)-FACT(3 ) ) =»STORA ( NNN , 1 
C0)/D£N 

STCRB(NE,11)=ST0RA(NNN, 11 )+NCRIT 

NB=NB+1 

GC TO 63 

61 ILDS=STCRA(NNN.ll )+NCRIT 

SPP.OB=FACT (2)*< 1 .0-FACT(l)-FACT(3) )*ST0RA(NNN,10) /DEN 
MCVES(ILOS)=MOVES( ILDS)+1 
PROB( ILCS)=PROB( ILCS)+SPROB 
63 DC 65 NNM=1,9 

MISSIUNNM)=STORA(NNN,NNM) 

65 CCNTINUE 
MISSIL(3J=0 
NCRIT=0 

CALL RELOAD ( M I SS I L , NCRIT ) 

IF (MISSIL(9).EQ.O) GO TO 69 

DC 67 IN=1,9 

STCRB(NB, IM=MISSIL(IN) 

67 CCNTINUE 

STORB(MB,10) = FACT(3)=4=(1.0-FACT( 1)-FACT( 2) ) *STORA( NNN , 1 
C0)/DEN 

STCRB(NB, 11)=ST0RA(NNN, 11 )+NCRIT 

NB=NB+1 

GO TO 70 

69 ILCS=STCRA(NNN,11)+NCRIT 

SFRCB=FACT(3)*(1 .0-FACT(l )-FACT( 2) ) =!‘STORA ( NNN , 10 ) /DEN 
MCVES(ILDS)=MOVES( ILDS)+1 
PRCB(ILCS)=SPR08+PPGE(ILDS) 

70 DC 71 LI=1,9 
M1SSIL(L1 j=STORA(NNN,LI) 
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\ 3)^^ ... ^^S i' i|j^= . *'l 1 : ' 

'■ " '"' " II '■ ’* " -Im '^'f 







71 CCNTINUe 
MISSIL(2)=0 
f<ISSIL(3)=0 
NCRIT=0 

C^LL RELOAD (MISS I L ,NCRIT ) 

IF (MISSIL(9) .EQ.O) GO TO 75 
DC 73 LN=1,9 
STORB(NE,LN)=MISSIL{LN) 

72 CCNTINUE 

STOPS (NB, 10)=2.0=i'FACT(2)*FACT(3)«(1.0-FACT( 1) )*STOPA(N 
CNN,10)/CEN 

STCRB(NB, ll) = STORA(NNNf ID+NCRIT 

NB=NB+1 

GC TO 77 

75 ILDS=STCRA(NNN, 11 I+NCRIT 

SPR0B=2 .0#FACT( 2 ) *FACT (3 ) * ( 1 .O^FACT( 1 ) ) 4=STCFA (NNN , 10 ) / 
CCEN 

MCVES( ILDS )=MOVES( ILDS )+l 
PROS (I LCS)=SPROB+PPCB (ILDS) 

77 DC 79 L2=l,9 

MSSILI LZ)=STCRA(NNN,LZ) 

79 CCNTINUE 
MISSIL( l)-0 
MISSIL( 3)=0 
NCR1T=0 

CALL RELOAD ( M IS S I L , NCRIT ) 

IF (MISSIL(9) .EQ.Ol GO TO 83 
DC 81 LA=1,9 
STCPB(N6 ,LA)=MISSIl(LA) 

81 CCNTINUE 

STCR3(Ne, 10)=2.0=*'FACT( 1)=J'FACT(3)*(1.0-FACT( 2) )*STCRA(N 
CNNf 10) /DEN 

STCRB(N8,11)=ST0RA(NNN,11)+NCRIT 

NB=NB+1 

GC TO 35 

82 ILCS=STCRA(NNN, 11 )+NCRIT 

SPRCB=2.0«FACT(1 )«FACT(3)=«'( l.C-FACTIZ) )*STCPA(NNN,10)/ 
CCEN 

«CVES(ILDS)=MOVES( ILDS)+1 
PR0B(ILCS)=SPP03+PPCe (ILDS) 

85 DC 87 LB=1,9 

MISSIL(L6)=ST0RA(NNN,LB) 

87 CONTINUE 
MISSIL( 1)=0 
MISSIL(2)=0 
NCRIT=0 

CALL RELOAD ( M ISS I L , NCR IT ) 

TO 91 

) 

1)*FACT(2 )*(1.0-FACT(2) )*STCRA(N 
,11)+NCRIT 



(2)«(1.0-FACT(3) )*STOPA(NNN, 10)/ 

) + l 

SPROB 

0 ) 



RIT) 

TO 99 

) 



IF (MISSIL(9) .EQ.O) GO 
CC 89 LC=1,9 
STCRe(Ne,LC)=MISSIL(LC 
89 CCNTINUE 

STCRB(NBf 10)=2.0=i‘FACT( 
CNN, 10) /DEN 
STCRB(NB,11)=ST0RA(NNN 
NB=NB+1 
GC TO 93 

91 ILCS=STCRA(NNN,11)+NCR 
SPP0B=2 .O^FACTd ) *FACT 
CCEN 

MCVES( ILDS )=MOVES ( ILDS 
PROB(ILDS)=PROB( I LDS)+ 
93 CC 95 LD=1,9 

MISSIL(LDJ=STORA(NNN,L 
95 CONTINUE 
MISSIL( I )=0 
MISSIL(2)=0 
MISSIL( 3)=0 
NCRIT=0 

CALL RELOAD (MISSIL,NC 
IF (MISSIL(9).EQ.O) GO 
CC 97 LE=1,9 
STCRB(NB,LE)=MISSIL (LE 
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c 

c 

c 

c 



S7 CONTINUE 

STCRBINE, 10) = 6.0*FACT( 1)*FACT(2)*FACT(3)=«'STCRA(NNN,10) 
C/DEN 

STCRB(NB,ll)=STORA(NNN, ID+NCRIT 

KB*NB+1 

GC TO 5000 

99 ILDS=STCRA(NNN»ll )+NCRIT 

SFRCB=6 .O^FACTIl) *FACT(2)*FACT(3)*ST0RA(NNN,10)/0EN 
MOVES! ILOS)=MGVES( IL0S)+1 
PRCe(ILCS)=PR03( ILDS)+SPROB 
GC TO 5C00 

TFE FORMULAE IN TABLE II ARE ENCODED BETWEEN 
HERE AND LINE 270. 

150 DEN=FACT(1)+FACT(2) 

MISSIL(2)=0 

NCRIT=0 

CALL RELOAD (MISSIL,NCRIT) 

IF (MISSIL(9J.EQ.O) GO TO 154 
DC 152 NA=1»9 
STCRB(NB,NA)=MISSIL(NA) 

152 CONTINUE 

STCRB(NB,10)=FACT (2)*( 1.0-FACT (1) ) *STORA ( NNN , 10 ) /DEN 

STORB(Ne,.ll)=STORA(NNN,ll)+NCRIT 

NB*NB+1 

GC TO 156 

154 ILDS=STCRA(NNN,11)+NCRIT 

SPR 0 e=FACT( 2 )*( 1 . 0 -FACT(l) ) =i'STC RA (NNN , 10 ) /DEN 
MCVES(ILDS)=MOVES( ILDS)+1 
PROB(ILCS)=PROB( ILCS)+SPROB 
156 DC 158 KB=lt9 

MISSI L( KB)»STCRA(NNN,KB) 

158 CONTINUE 

MISSIL(2)=0 

NCRIT=0 

CALL RELOAD ( M ISS I L , NCR IT ) 

IF (MISSIL(9) .EQ.O ) GO TO 162 
DC 160 KC=lt9 
STCRS<NE,KC)=MISSIL(KC) 

16C CONTINUE 

STCRB(NB,10) = FACT( 1)*(1.0-FACT( 1) -FACT! 2) )=»STCRA! NNN, 1 
COJ/DEN 

STCRB!MB,11)=ST0RA!NNN,11)+NCRIT 

NB=NB+1 

GO TO 164 

162 ILDS=ST0RA!NNN.11)+NCRIT 

SFROB = FACT! 1)=!'! 1.0-FACT!!) -FACT! 2) )*STORA!NNN, 10) /DEN 
MCVES!ILDS)=MOVES! IL0S)+1 
PROB!ILCS)=SPROB+PROB!ILDS) 

164 DC 166 K0=lt9 

MISSIL!KD)=STORA!NNN,KD) 

166 CONTINUE 

MISSIL! 1)=0 
MISSIL! 3)=0 
NC.=!IT=0 

CALL RELOAD ! MI S S I L , NCR IT ) 

IF !MISSIL!9).EQ.O ) GO TO 170 
DC 168 KE=1»9 
STCRB!N6,KE)=MISSIL!KE) 

168 CONTINUE 

STCRB!NE,1C)=FACT!1 )**2*!1. 0-FACT !2) )*STGRA !NNN,10)/CE 
CN 

STrRB!NE,ll)=STORA!NNN, 11 )+NCRIT 

NE=NB+1 

GC TO 172 

170 ILDS=STCPA!NNN,11)+NCRIT 

SPPCB=FACT!l)*=i<2«! 1.0-FACT ! 2) )*STCRA!NNN, 10) /DEN 
MOVES! ILDS) = MOVES( ILDS)+1 
PROB!ILDS)=PPCB( 1LDS)+SPR0B 
172 DC 174 KF=lt9 

MISSIL! KF) =STCRA!NNN,KF) 
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174 CCNTINUE 

«ISSIL(2)=0 

MISSIL(3)=0 

NCRIT=0 

CALL RELOAD (MISSIL,NCRITJ 
If (MISSIL(9) .EQ.O ) GO TO 178 
DC 176 KG=1,9 
STCR8(NBtKG)=MISSIL(KG) 

176 CCKTINUE 

STCRB(NB, 10)=2.0*FACT(1)«FACT(2)’«‘(1.0-FACT( in*STORA(N 
CNN*10)/DEN 

STCR6(NBtll)=ST0RA(NNNtll)+NCRIT 

N6=NB-H 

GC TO 180 

176 ILDS=STCRA(NNN,ll)-«-NCRIT 

SPR08=2.0=i'FACT( 1)«FACT(2)*( 1.0-FACT( in=!‘STCPA(NNN,10)/ 
CDEN 

MCVES( ILDS)=MCVES( ILDS)+1 
PROB(ILCSJ=PROB( ILCS)-«-SPROB 
180 DC 182 KH=1,9 

RISSIL(KH)=STGRA(NNN,KH) 

182 CONTINUE 

MISSIL( 1)=0 
MI$SIL( 2)=0 
MISSIL(3)=0 
NCRIT=0 

CALL RtlOAD (MISS IL ,NCRIT ) 

IF (M1SSIL(9) .EQ.O) GO TO 186 

DO 184 KI=lt9 I 

STGRB(NE,KI)=MISSIL(KI J 
184 CONTINUE 

STCRB(NBt 10)=3.0«FACT( 1)=!'*2*FACT(2)=5'ST0RA(N'NN,10) /DEN 
STCR6(Ne,ll)=ST0RA(l\NN, 11) + NCRIT 
NE=NB-*-l 
GC TO 5000 

186 IL0S=STCRA(NNN,11 )+NCRIT 

SPR0B=3.0=!'FACT( l)*=S'2-FACT(2)*STCBA(NNNt 10)/CEN 
MCVES( ILCS)=M0VES( ILDS)+1 
PRCBd LCS)=PROB(ILCS) + SPROB 
GO TO 5000 
C 

C THE FORMULAE IN TABLE III ARE ENCCDED BETbEEN 

C HERE AND LINE 390. 

C 

270 DEN=FACT(1)+FACT(2) 

MISSIL( 1)=0 
NCRIT=0 

CALL RELOAD ( MI SS I L ,NCRIT ) 

IF (MISS1L(9) .EQ.O) GO TO 274 
DC 272 IA=lt9 
STCR6(Nfi,IA)=MISSIL( lA) 

272 CCKTINUE 

STCRB(NB,10)=FACT( 1)*(1.0-FACT(2) )*ST0RA(KKh,10)/DEN 
STCRB(NE,11)=ST0RA( KNN, 11)+NCRIT 
NB=NE-H 
GO TO 276 

274 ILDS=STCRA(NNN,11)+NCRIT 

S PROB=F ACT ( 1 ) * ( 1 . 0-F 6CT ( 2 ) ) *STOR A ( NNN , 1 0 ) /D EN 
MCVES( ILDS )=MOVES( ILDS)+1 
FBC8(ILCS)=PR06( ILCS)+SPROB 
276 DC 27S IB=1,9 

MISSILC IB)=STCRA(NNN,IB) 

278 CONTINUE 

MISSIL(3)=0 

KCRIT=0 

CALL RELOAD ( M I S S I L , NCR IT ) 

IF (MISSIL(9) .EQ.O ) GO TO 282 
DC 280 IC=1,9 
ST0RB(N6,IC)=MISSIL( IC) 

280 CCNTINUE 

STCRBI NBtlO)=FACT (2)*(1.0-FACT( 1)-FACT(2) ) *STCRA( NNN » 1 
CO)/DEN 
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STCRB(NB,11)=ST0RA(NNN,11)+NCRIT 

NB=NB+1 

GC TO 2 84 

282 ILDS=STGRA(NNN,11 )+NCRIT 

SFBOB=FACT(2)«( 1 . 0-FACT ( 1 ) -FACT ( 2 )) *STORA ( NNN, 10 ) /DEN 
MCVES( ILOS)=MOVES( ILDS)+1 
PPOB(ILCS)=PROB( ILDS)+SPROB 
284 DC 286 IE=lt9 

NISSIL( IE)=STGRA(NNN,IE) 

286 CCwTINUE 

MISSIL(1)=C 

MISSIL(3)=0 

NCRIT=0 

CALL RELOAD ( M IS SI L ,NCRIT ) 

IF (MISSIL(9).EQ.O) GO TO 290 
DC 288 IF=1,9 
STCRB(NB,IF)=MISSIL(IF) 

288 CONTINUE 

STCRB (NE, 101=2. 0*FACT(1)=!'FACT (2 )«( 1.0-FACT ( 2) )*STORA(N 
CNKj,10)/DEN 

STCRB(NE,I1)=ST0RA(NNN,11 I+NCRIT 

NE=NB+1 

GC TO 292 

290 ILCS=STCRA(NMM,11)+NCRIT 

SPR06=2.0*FACT( 1)*FACT(2)*( 1.0-FACT( 2) ) *STC PA ( NNN , 10 ) / 
C C E ^ 

MCVES(ILDS)=NOVES( ILDS)+1 
PRCB(ILCSJ=PRCB( ILDS)+SPR06 
292 CC 294 IG=1»9 

RISSIL( IG)=STORA(NNN,IG) 

294 CONTINUE 

MISSIL(2)=0 

MISSIL(3)=0 

NCRIT=0 

CALL RELOAD ( M I S S I L ,NCR IT ) 

IF (MISSIL(9) .EQ.O ) GO TO 298 
DC 296 IH=1,9 
STCRB(Nfi,Ih)=MISSIL(IH) 

296 CCKTINUE 

STCRB(NE,10)=FACT(2)**2*(1.0-FACT(1))=«=STCRA(NNN,10)/CE 

CN 

STCRB(NE,11)=STCRA(NNN, 11 J+NCRIT 

N6=NB+1 

GC TO 300 

298 ILDS=ST0RA(NNN,11)+NCRIT 

SPROB=FACT(2)**2#( l.C-FACT( 1) )*STCRA(NNN , lO/DEN 
MOVES ( ILDS )=MOVES( ILDS )+l 
PRCB( ILCS)=PROE( ILDS)+SPROB 
300 DC 302 IQ=lt9 

MISSIL( IQ)=STORA(NNN,IQ) 

302 CONTINUE 

KISSIL(1)=0 
MISSIL( 2)=0 
MISSIL(3)=0 
NCPIT=0 

CALL RELOAD ( M I SS I L , NCR IT ) 

IF (MISSIL(9).EQ.0 ) GO TO 306 
DC 304 IR=1,9 
STCRB(NB,IR)=MISSIL(IR) 

304 CONTINUE 

STCRB(NB,10)=3.0=S'FACT(1)=!=FACT(2)**2*STGRA(NNN,10)/DEN 
STCRB(NB,11)=ST0RA(NNN, ID + NCRIT 
NE = N'B+1 
GC TO 5000 

306 ILCS=STCRA(NNN,11)+NCRIT 

SPR0B=3.0*FACT(1 )«FACT(2)**2*ST0PA(NNN, 10) /CEN 
MCVES(ILDS)=MCVES( ILDS)+1 
FRCB(1LCS)=PR0B( ILCS)+SPROB 
GC TO 5000 
C 

C THE FORMULAE IN TABLE IV ARE ENCCCED BETWEEN 

C HERE AND LINE 500. 
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c 

390 DEN=FACT(1)+FACT(3) 

MISSIL(2)=0 

NCRIT=0 

CALL RELOAD ( M I SS I L , NC R IT ) 

IF (MISSIL(9).EQ.O) GO TO 394 
00 392 ^'A = 1,9 
STCRB(Ne,MA)=MISSIL(MA) 

392 CONTINUE 

STGRB(NB»10i=FACT(l)«(1.0-FACT( 1)-FACT(3) ) *STCRA ( NNN , 1 
CO)/DEN 

STCRB(NB,I1)=ST0RA(NNN,11 )+NCRIT 

NB=NB+1 

GC TO 396 

394 ILDS=ST0RA(NNN,11)+NCRIT 

SFROB=FACT( 1)*( 1.0-FACT( 1)-FACT( 3) )-ST0RA(NNN,10) /DEN 
MCVES( ILDS)=MOVES( ILDS )+l 
PBCB(ILCS)=PRGB( ILDS)+SPROB 
396 DO 398 IC=li9 

PISSIL( IC)=STCRA(NNN',IC) . 

398 CONTINUE 

PISSIL(3)=0 

NCRIT=0 

CALL RELOAD ( MI SS IL , NCR IT ) 

IF (MISSIL(9).EQ.O) GO TO 402 
CO 400 ID=1,9 
STCRB(NBfID)=MISSlL( ID) 

400 CONTINUE 

ST0RB(NB,10)=FACT(3)=i=( 1.0-FACT( 1) ) *STQR A ( NNN , 10) /DEN 

STCRB(NBtll)=STORA(NNN,ll)+NCRIT 

KB=NB+1 

GC TO 404 

402 I LDS=STCRA(NNN,11 )+NCRIT 

SPR0B=FACT(3)*a .C-FACT( 1) )«ST0PA(N'NN»10) /DEN 
MCVES( ILDS)=M0VES( ILDS)+1 
PROBd LCS)=PROB( ILCS)+SPRQB 
404 DC 406 IE=1,9 

MSSIL( IE)=STORA(NNN,IE) 

406 CONTINUE 

MISSIL( 1)=0 
MISSIL(2)=0 
NCRIT=0 

CALL RELOAD (MISS IL ,NCRIT ) 

IF (MISSIL(9) .EQ.O) GO TO 410 
CC 408 IF=1,9 
STGRB(NB,IF)=MISSIL( IF) 

408 CONTINUE 

ST0RB(NE,10) = FACT( 1 )**2*(1.0-FACT( 3) )=4'STCPA(NN'N,10)/CE 
CN 

STCPB(NB»11)=ST0PA(NNN,11) + N'CRIT 

NB=NB+1 

GO TO 412 

410 ILCS=ST0RA(NNN,11 )+NCRIT 

SFP0B=FACT(l)=^•*2=i*( 1.0-FACT( 3) ) ’i^STC PA ( NNN , 10 )/DEN 
MOVES ( I LDS)=MGVES( IL0S)+1 
PPOB(ILCS)=PROB( ILDSJ+SPROB 
412 CC 414 IG=1,9 

MISSIL( IG)=STORA(NNN,IG) 

414 CONTINUE 

MISSIL(2)=0 
MISSIL( 3)=0 
N'CRIT=0 

CALL RELOAD ( M I SSI L ,NCR IT ) 

IF (MISSIL(9) .EQ.O) GO TO 418 
DC 416 IH=1,9 
STCRB(Ne,IH)=MISSIL(IH) 

416 CONTINUE 

S7CRB(N6, 10) =2.0*FACT( 1 ) *F ACT ( 3 ) =!= ( 1 .0-FACT ( 1 ) )*STOPA (N 
CNN,10)/DEN 

STCPB(NE,11)=ST0RA(NNN,11)+NCRIT 

NB*NB+1 

GC TO 420 
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418 ILDS=STCR-A(NNN,11)+NCRIT 

SPRCB=2.0*FACT(l)*FACT(3)«(1.0-F/iCT(l ))*STCPA(NNN,10)/ 
COEN 

MCVES( ILOS)=MOVES( ILDS)+1 
PBOB(ILOS)=PROB(ILDS)+SPROB 
420 DC 422 IJ=1,9 

RISSIK IJ)=ST0RA(NNN, IJ) 

422 CONTINUE 

KISSIU 1)=0 
MISSIL(2)=0 
MISSIL(3)=0 
NCPIT=0 

CALL RELOAD ( MI S S I L ,N CRIT ) 

If (MISSIL(9) .EQ.O J GO TO 426 
DC 424 IK=1,9 
STGRBINB, IK)=MISSIL(IK) 

424 CONTINUE 

STCRB(NB,10)=3.0*FACT(1)**2*FACT(3)*STORA(NNN,10)/DEN 

STCRB(NB,11}=STQRA(NNN,11)+NCRIT 

NB=NB+1 

GO TO 5 COD 

426 IL0S=ST0RA(NNN,11)+NCRIT 

SFRGB=3.0^FACT( 1)=!'*2=!'FACT(3)*STGRA(NNN,10)/CEN 
MOVES! ILOS)=MOVES( ILDS J+1 
PPCB(ILDS)=PROB( ILCSJ+SPROB 
GO TO 5000 
C 

C TFE FORMULAE IN TABLE V ARE ENCODED BETWEEN 

C HEPE AND LINE 5000. 

C 

500 MISSIL(3)=0 
NCPIT=0 

CALL RELOAD ( MI S S I L ,NCR IT ) 

IF (MISSIL(9) .EQ.O) GO TO 504 
CG 502 JQ=lt9 
STCRBCNE, JQ)=MISSIL( JQJ 
502 CONTINUE 

STGRB(N6,10)=( l.O-FACT! IJ ) *STOR A ( NNN , 10 ) 
STGRB(NB,11)=ST0RA(NNN,11 )+NCRIT 
NB=NB+1 
GO TO 506 

504 ILDS=STORA(NNNtll J+NCRIT 

SFPOB=( 1.0-FACT( 1) ) *STCRA ( NNN , 1 0 * 

MGVES(ILDS)=MOV£S( ILDSJ+1 
PPOB (I LOS )=PROe( ILCSJ+SPROB 
506 00 508 JR=1,9 

MISSIL! JR) =STCRA(NNN,JR) 

508 CONTINUE 

MISSIL! 2)=0 
MISSIL(3)=0 
NCPIT=0 

CALL RELOAD ( MIS S I L , NC R IT J 
IF (MISSIL(9) .EQ.O) GO TO 512 
00 510 J5=lt9 
STCPBINB, JS)=MISSIL( JS) 

510 CONTINUE 

ST0RB(NB,10)=FACT( 1J=^( 1.0-FACT{ 1) ) 4ST0RA ( NNN , 10 ) 

ST0RB(NE,I1)=STQRA(NNN,11)+NCRIT 

NB=NB+1 

GO TO 514 

512 ILDS=STQRA(NNN, 11 J+NCRIT 

SPRCB=FACT(1)*( l.O-FACT(l) ) *STOR A ( NNN , 10 ) 

MOVES! ILDS)=MOVES! ILOS)+i 
PRCB ! I LDSJ=PROB! ILCSJ+SPROB 
514 DC 516 JT=i,9 

MISSIL! JTJ=STORA!NNN,JTJ 
516 CONTINUE 

MISSIL! 1 J=0 
MI5SIL! 2J=C 
MISSIL!3J=0 
NCRIT=0 

CALL RELOAD ! MI S S I L ,NCR IT J 
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IF (MISSIL(9) .EQ.O) GO TO 520 
DC 518 JU=1,9 
STCRB(Ne,.JU)=MISSIL( JU) 

518 CONTINUE 

STCRB(NE,10)=FACT( 1 )«*2*ST0RA ( NNN , 10 ) 
STGRB(NB,11)=ST0RA (NNN, 11 )+NCRIT 

Na=N8+l 

GC TO 5000 

520 ILDS=STCRA(NNN,11 I+NCRIT 

SPRCB=FACT(1)*=<'2=*'STGRA(NNN, 10) 

MCVESI ILDS)=M0V£S( ILOS )+l 
PRGBI ILCS)=PR08( ILDS)+SPROB 
GC TO 500.0 
5000 CONTINUE 

THIS OUTPUT GIVES THE MASS FUNCTION OF THE NUMBER 
OF MOVES REQUIRED TO EMPTY THE STORAGE BARNS AND 
ALSO THE NUMBER OF BRANCHES YIELDING EACH PARTICULAR 
VALUE CF NUMBER OF MOVES. 

6000 WRITE (6,8095) 

8095 FORMAT ('O', 'THE PROBABILITIES OF MOVES ARE THE OCCUR 
CRENCES OF MOVES ARE' ) 

DO 600 1=1,40 

WRITE (6,8100) I,PPOB(I),MOVES{I) 

8100 FORMAT (• • , 13 , 5X , FIO . 8 , lOX , 1 10 ) 

600 CONTINUE 
IPATH=0 
PSUM=0.0 
DO 700 JJ=1,40 
PSUM=PSUM+PR0B( JJ ) 

IPATH=I PATH+MOVES ( J J ) 

7C0 CONTINUE 

THIS OUTPUT ENABLES A VISUAL CHECK TO DETERMINE IF 
THE MASS FUNCTION SUMS TO ONE AND HOW MANY BRANCHES 
EXISTED IN THE ENTIRE TREE. 

WRITE (6,8105) PSUM,IPATH 

8105 FORMAT ('O', 'SUM CF PROBABILITIES I S • , FIO . 8 , • TOTAL PAT 
CHS COUNTED* ,110) 

STOP 

END 



SUBROUTINE RELOAO( LNCHRS, MOVES ) 

DIMENSION LNCHRS(9) 

THIS SUBROUTINE RELOADS THE EMPTY LAUNCHERS AND COUNTS 
THE NUMBER OF MOVES REQUIRED. WHEN THE BARN IS EMPTY, 
THE FINAL TOTAL OF MOVES IS OUTPUT. THE RELOAD 
DOCTRINE ilS TO MOVE ONLY ENOUGH MISSILES TO REFILL 
THE LAUNCHERS. 

LCHEPT=0 
LST0RE=0 
DC 5 MI=1,3 

IF (LNCHRS(MI).EQ.O) LCHEPT=LCHEPT+1 
5 CONTINUE 
DC 10 MJ=4,9 

IF (LNCHRS(MJ) .NE.O) LSTORE=LSTCR E+1 
10 CONTINUE 

IF (LCHEPT. EQ.O. OR. LSTORE. EQ.O) RETURN 
IF (LCHEPT. LE. LSTORE) GO TO 70 
60 TO (15, 20, 40), LCHEPT 
15 RETURN 

2C IF (LNCHRSO) .NE.O) GO TO 33 
DC 25 MK=4,9 
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IF (LNCHRS(MK).NE.O) GO TO 30 
25 CONTINUE 

30 LNCHRS( 3)=LNCHRS(MK) 
LNCHRS(NK)=0 
FCVES=MCVES +(MK-3) 

RETURN 

33 LNCHRS(2)=LNCHRS(3) 

DC 35 MKK=4»g 

IF (LNCHRS(MKK) .NE.O) GO TO 37 
35 CONTINUE 

37 LNCHRS (3 ) =LNCHRS ( MKK ) 
LNCHRS(MKK)=0 
MOVES =MCVES+l+(MKK-3) 

RETURN 

^0 IF (LSTCRE.NE.l) GO TO 55 
A3 DC 45 ML=4»9 

IF (LNCHRS(ML).NE.O) GO TO 50 
45 CONTINUE 

50 LNCHRS(3)=LNCHRS(ML) 

LNCHRS ( ML )=0 
M0VES=M0VES+(ML-3) 

RETURN 

55 DC 60 MM=4»9 

IF (LNCHRS(MM).NE-O) GO TO 65 
60 CONTINUE 

65 LNCHRS( 2)=LNCHRS(MM) 

LNCHRS( MM) =0 
MCVES=M0VES+(MM-2) 

60 TO 43 
70 DO 100 MN=1,3 

IF (LNCHRS(MN).NE.O) GO TO 100 

MIN=MN+1 

DC 75 MII=MINt9 

IF (LNCHRS(MII) .NE.O) GO TQ 80 
75 CONTINUE 

80 LNCHRS(MN)=LNCHRS(MII ) 
LNCHRSCMLI ) = 0 
MOVES=MCVES+(MII>MN) 

100 CONTINUE 
RETURN 
END 
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COMPUTER PROGRAM II 



DIMENSION MISSIL(9J ,PR0P0R(5), INARAY(9) ,KTGTCT(5) , 
CSLIMIT(A)»PWRHDS(5) , P MOVE ( 40) , J MOVES ( 40) 

DATA KSEED/l91/iKTGTCT/5'-!'0/,JM0VES/40*0/, WRhDS/0.0/, 
CLGCKON/O/,RMOVES/0.O/ 

READ IN THE NUMBER CF MISSILES I^ THE SECTION. 

READ (5,1050) IBASLD 
1050 FORMAT (• *,11) 

READ IN THE INITIAL ORDER OF THE WARHEADS. 

READ (5,1200) INARAY 
1200 FORMAT (912) 

READ IN THE PROBABILITIES OF A REQUEST FOR A 
PARTICULAR WARHEAD. THE FIRST PROBABILITY IS FOR THE 
TYPE 1 WARHEAD, THE SECOND PROBABILITY IS FOR THE TYPE 
2 WARHEAD,.ETC. FIVE NUMBERS MUST BE INPUT. . THE FIVE 
NUMBERS MUST ADD TO ONE. 

READ (5,1300) PROPCR 
1300 FORMAT (5F10.3) 

THIS OUTPUT ENSURES THAT THE INITIAL ARRAY WAS INPUT 
CORRECTLY. 

WRITE (6,1400) INARAY 

1400 FORMAT (• ','THE INITIAL ARRAY IS', 912) 

SLIM.IT( 1) = PR0P0R( 1) 

DC 10 J=2,4 

SLIMIT( J)=SLIMIT( J-1)+PR0P0R( J) 

10 CONTINUE 

THIS OUTPUT ENSURES THAT THE WARHEAD PROBABILITIES 
WERE INPUT PROPERLY. 

WRITE (6,1600) ( PROPOR(NNN) ,NNN=1,5) 

1600 FORMAT (• ','THE WARHEAD PROBABILITIES ARE',5F10.3) 

THIS SIMULATION WILL EE RUN 10,000 TIMES. EACH RUN 
WILL BE WITH A DIFFERENT RANDOM STRING OF REQUESTS 
FOR WARHEADS. 

DO 50 K = l,.10000 
DC 14 KK=1,9 
MISSIL(KK)=INARAY(KK) 

14 CONTINUE 
NCRIT=0 
NCCUNT=0 

15 CALL TARGET(SLIMIT ,LOCKON,KTGTCT,KSEED) 

DO 20 N=l,3 

THIS CHECKS IF A LAUNCHER CONTAINS A WARHEAD OF THE 
TYPE REQUESTED. THE LAUNCHERS ARE CHECKED STARTING 
WITH THE ONE CLOSEST TO THE STORAGE BARNS, THEN 
CHECKING THE NEXT LAUNCHER OVER, AND FINALLY 
CHECKING THE LAUNCHER FARTHEST FROM THE BARN. 



IF (MISSIL(4-N).EQ.L0CK0N) GO TO 30 
20 CONTINUE 

CALL RELOAD ( MIS S IL , NCRIT ) 

GO TO 15 
30 MISSIL( 4-N)=0 
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NCCUNT=NC0UNT+1 

IF (NCOUNT.lt. IBASLD) GO TO 15 
JN0VES(NCRIT)=JM0VES(NCRIT)+1 
50 CONTINUE 

DC 60 11=1,5 
WRHDS=WRHDS+KTGTCT (II) 

60 CONTINUE 

THIS OUTPUT ENABLES A VISUAL CHECK OF WHETHER THE 
RANDOM ST-RING OF REQUESTS CONTAINS THE PPCPER 
PROPORTION OF TYPES. 

WRITE (6,1800) 

1800 FORMAT (• •, ‘WARHEAD DISTRIBUTION IS*) 

DC 70 IJ=1,5 

PWRHOS( I J)=KTGTCT( I J)/WRH0S 
WRITE (6,1900) I J, PWRHDS( I J) 

1900 FORMAT (• • , 1 10 , 5X , FIO . 8 ) 

70 CONTINUE 

DC 75 IL=1,40 
RMOVES=RMOVES-fJMOVES( ID 
75 CONTINUE 

THIS OUTPUT PROVIDES THE MASS FUNCTION OF THE NUMBER 
OF MOVES ^^EOUIRED TO EMPTY THE STORAGE BARNS DURING 
THE 10,000 TRIALS. 

WRITE (6,2000) 

2000 FORMAT (• • , • DI STR I BUT I ON OF MOVES IS*) 

DO 80 IK=1,40 

PMOVE( IK)=JMOVES( IK)/RMOVES 
WRITE (6,1900) IK,PMOVE(IK) 

80 CONTINUE 
STOP 
END 



SUBROUTINE TARGET ( POINTS, MEXP, JCCUNT,IX) 

DIMENSION POINTS (4) ,JC0UNT(5) 

THIS SUBROUTINE GENERATES THE RANDOM STRING OF 
WARHEAD RtGUESTS. A RANDOM NUMBER BETWEEN ZERO AND 
0^E IS GENERATED AND THEN THIS NUMBER IS COMPARED 
WITH THE INPUT WARHEAD PROBABILITIES IN ORDER TO 
DETERMINE THE TYPE OF REQUEST GENERATED. 

IX=IX*65539 

YFL=0.5-HX*2.328306E-10 
DO 5 M=l,4 

IF (YFL.LE.PGINTS(M) ) GO TO 10 
5 CONTINUE 
MEXP=5 
60 TO 20 
1C MEXP=M 

20 JCCUNT(MEXP)=JC0UNT(MEXP)+1 
RETURN 
END 



SUBROUTINE RELOAD ( LNCHRS, MO VES ) 

DIMENSICN LNCHRS(9) 

THIS SUBROUTINE RELOADS THE EMPTY LAUNCHERS AND CCLNTS 
THE NUMBER OF MOVES REQUIRED. WHEN THE EARN IS EMPTY, 
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THE FINAL TOTAL OF MOVES IS OUTPUT. THE RELOAD 
DOCTRINE IS TO MOVE ONLY ENOUGH MISSILES TO REFILL 
THE LAUNCHERS. 

LCHEPT=0 
LST0RE=0 
DO 5 MI=1,3 

IF (LNCHRS(MI ) .EQ.O) LCHEPT=LCHEPT+1 
5 CONTINUE 
DC 10 MJ=4,9 

IF (LNCHRS(MJ) .NE.O) LSTORE=LSTOR E+l 
10 CONTINUE 

IF (LCHEPT. EQ.O. OR. LSTORE .EQ.O) RETURN 
IF (LCHEPT. L£. LSTORE) GO TO 70 
GC TO ( 15,20,40) , LCHEPT 
15 RETURN 

20 IF (LNCHRSO) .NE.O) GO TO 33 
DC 25 MK=4,9 

IF (LNCHRS(MK).NE.O) GO TO 30 
25 CONTINUE 

30 LNCHRS(3)=LNCHRS(MK) 

LNCHRS(MK)=0 
MCVES=MOVES -KMK-3) 

RETURN 

33 LNChRS( 2)=LNCHRS(3) 

DC 35 MKK=4,9 

IF (LNCHRS(MKK) .NE.O) GO TO 37 
35 CONTINUE 

37 LNCHRS (3 ) =LNCHRS ( MKK ) 

LNCHRS( MKK)=0 

MOVES =MQVES+l+(MKK-3) 

RETURN 

40 IF (LSTORE. NE.l) GO TO 55 
42 DC 45 ML=4,9 

IF (LNCHRS(ML) .NE.O) GO TO 50 
45 CONTINUE 

50 LNCHRS(3) = LNCHRS( ML) 

LNCHRS( ML) -0 
MCV£S=MCV£S+(ML-3) 

RETURN 

55 DC 60 MiM=4,9 

IF (LNCHRS(MM) .NE.O) GO TO 65 
60 CONTINUE 

65 LNCHRSI 2)=LNCHRS(MM) 

LNCHRS ( MM) =0 
MGVeS=M0VES+(MM-2) 

60 TO 43 

70 DC 100 MN=1,3 

IF (LNCHRS(MN) .NE.O) GO TO 100 

MIN=MN+1 

DC 75 MII=MIN,9 

IF (LNChRS(MII) .NE.O) GO TO 80 
75 CONTINUE 

80 LNCHRS( MN)=LNCHRS(MII ) 

LNCHRS(MM ) = 0 
MCVES=MCVES+(MI I-MN) 

100 CONTINUE 
RETURN 
END 
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